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Within the imaginary-time theory for nonequihbrium in quantum dot systems the calculation of 
dynamical quantities hke Green's functions is possible via a suitable quantum Monte-Carlo algo- 
rithm. The challenging task is to analytically continue the imaginary-time data for both complex 
voltage and complex frequency onto the real variables. To this end a function-theoretical descrip- 
tion of dynamical observables is introduced and discussed within the framework of the mathematical 
theory of several complex variables. We construct a feasible maximum-entropy algorithm for the 
analytical continuation by imposing a continuity assumption on the analytic structure and provide 
results for spectral functions in stationary non-equilibrium and current-voltage characteristics for 
different values of the dot charging energy. 

PACS numbers: 



I. INTRODUCTION 

Dynamic observables play a crucial role in 
the Matsubara voltage approach introduced by 
Han and Heary to address steady-state non- 
equilibrium properties of models for quantum 
dots.fil In a previous publication (hereafter re- 
ferred to as Ijpl we showed under what con- 
ditions the real-time Keldysh and imaginary- 
time Matsubara-voltage approaches are formally 
equivalent, and how a proper analytical contin- 
uation must be performed. In I, this scheme 
was applied to static quantities obtained from 
a continuous-ti me q uantum Monte-Carlo (CT- 
QMC) algorithnplifil combined with a standard 
maximum-entropy (MaxEnt) approactP^ to ob- 
tain results for steady-state expectation values of 
quantum dot models at finite bias. 

We consider a single-impurity Anderson model 
for the quantum dot systerrP^ with the Hamilto- 
nian 

cxka <y 
-E4=(4co..+h.C.) (1) 
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This Hamiltonian describes a quantum dot device 
which consists of the quantum dot orbital opera- 
tor of spin a and with source and drain leads, 
represented by conduction electron operators cl^i^^ 
with the continuum index fc, the spin index a and 
the reservoir index a = L,R for the source (left) 
and drain (right), respectively. The model is char- 
acterized by few parameters;'*" The local energy or 



gate voltage Cd, which controls the number of elec- 
trons on the dot; the charging energy U > due 
to the small capacitance; and finally the coupling 
of the dot to the leads, which can in many cases 
be collected in two quantities, and F^j. 

Although some thermodynamic observables can 
be calculated directly without analytic continu- 
ation of imaginary-frequencj- ^ , only a restricted 
set of observables can be handled in this man- 
ner. Unfortunately, the current operator (/) is 
not suitable, and the rather important question 
about the transport through a quantum dot, both 
electrically and thermally driven, must be ad- 
dressed in a different manner. For simple quan- 
tum dot geometries one can employ the result by 
Mcir and Wingreen,'^ who showed that for single 
quantum dots and not too different properties in 
the left and right leads, one can express the cur- 
rent through the dot due to a finite external bias 
cVb = $ via the density-of-states (DOS) 7^4(6; $) 
on the dot as 

m) =hY.lde [fde) - fR[e)] M.[e- $) (2) 

where /«(£) denotes Fermi's function for the left 
respectively right lead. Note that the bias enters 
in two distinct ways: First, in the Fermi func- 
tions through the chemical potential of the leads 
afP^/x±$/2, and second, through the DOS. Usu- 
ally, for <I> — >■ 0, one ignores the latter dependency 
and can thereby recover the results from linear- 
response theory. 

Besides its relevance for calculating the current, 
the DOS is an interesting quantity in its own right, 
and its dependency as function of frequency and 
bias is still a matter of debate. In equilibrium, it 
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is well-known that the DOS develops a very sharp 
resonance, the so-called Kondo resonance, pinned 
at e = 0, which in linear response leads to the 
pinning of the conductance to the unitary limit. 
The precise way that this resonance dies under the 
influence of finite bias is actually unknown, and 
different techniques provide different answers. 

Within the Matsubara-voltage imaginary-time 
approach by Han and Heary the nonequilibrium 
steady state is mapped to an infinite set of effec- 
tive equilibrium systems by intro ducing a bosonic 
Matsubara voltage ip^ = 47rTO//3.'32lIIi2l It has to 
be analytically continued to a variable z,^ € C in 
order to compute the limit ^ $ ± W to obtain 
the physical quantity at the chosen bias $. 

Considering dynamic observables, the simulta- 
neous presence of the fermionic Matsubara fre- 
quency iaj„, which must be analytically contin- 
ued to g C to eventually obtain the corre- 
sponding real-frequency quantity, implies that a 
double-complex- variable function G(z), with z := 
[.z^p^z^y G C^, must be studied. Without the 
presence of it^m, i.e. within conventional Matsub- 
ara theory. Green's functions are analytic on the 
upper and lower half-planes, IH and H*. Due to 
the rapid decay of G(ztj) as z^j — > cx), a Lehmann 
spectral representation with respect to the real 
axis is used.I^l Conversely, from a knowledge of all 
Matsubara- frequency data one can in principle in- 
fer the spectral function. Numerically, this is a 
known to be an ill-conditioned problem. An ap- 
proach, particularly suited for QMC data, is the 
MaxEnt technique.!^ 

From a mathematical point of view, the branch 
cut on the real axis represents the only set of 
points zo G C for which the Green's function is 
not holomorphic. The very location of the branch 
cut gives rise to the spectral representation in 
Matsubara theory, i.e. it yields, due the nonsin- 
gular structure at oo, a unique characterization 
of the analytic structure. In order to perform an 
analytic continuation for the two-variable func- 
tion G(z) the natural question arises as to which 
minimal set of quantities characterizes its analytic 
structure in a unique fashion. Therefore the re- 
ferral to mathematical results in the analysis of 
several complex variables is necessary (see in p ar- 
ticular the discussion following section III C 1 ).^" 



In the original implementation of Han and his 



co-workepi^, analytic continuation on (z,^,^tj; 
produced rather smooth spectra in g ood agree- 
ment with other numerical result d^^^^l However, 
the employed methods were crude fit with ad- 
hoc smoothening and annealing to the numerical 
data without regards to statistical analysis. In 
this work we provide mathematical foundation to 



multi-variable analytic continuation and develop 
systematic numerical implementation. 

The paper is structured as follows. Since the 
mathematical structure dealing with functions of 
several complex variables is probably very alien 
to the reader, we start with a presentation of the 
central results for spectral functions and transport 
properties of the single-impurity Anderson model 
in steady-state non-equilibrium in section [Hj The 
underlying mathematical framework is developed 
in the succeeding sections. Starting from conven- 
tional one-dimensional complex analysis, section 
JH] provides an introduction to the basic concepts 
of theory of several complex variables. We then 
use this theory to systematically analyze the ana- 
lytical structure of the Matsubara-voltage Green's 
function and provide an axiomatic description of 
it. In section IV an asymptotically exact conti- 



nuity assumption is employed to construct an in- 
tegral representation for Matsubara Green's func- 
tion G(i(p,„, iw„). It allows us to include more in- 
formation within the Bayesian inference process of 
the maximum entropy method (MaxEnt) as com- 
pared to our previous approach in Ref. 1101 The 
resulting structure connects to the earlier sugges- 
tions by Han and Heary.^i^ For future applica- 
tions, an unbiased extension of the integral rep- 
resentation beyond the continuity assumption is 
proposed in section [V] 

Figure [l] provides an overview of the paper as 
a flowchart. It may serve as a guide. Mathemat- 
ically less inclined readers may skip the lightly 
shaded section [ffl] and study the numerical re- 
sults in section |TT| which are based on the con- 
tinuity ansatz and MaxEnt procedure described 
in section IV and corresponding appendices. In 
the chart full arrows denote numerical steps in 
the computation. Empty arrows denote formal 
analytical steps required to derive the relations 
involved into the MaxEnt procedure for analytic 
continuation. 

All numerical results provided in the paper 
rely on a highly precise continuous-time quantum 
Monte-Carlo (QMC) implementation which was 
introduced in Ref.[TOl The data provided by QMC 
simulations, see left box in figure [T] give rise to 
a discrete grid of well-estimated imaginary-time 
Green's function values. Due to the structure of 
the Matsubara-voltage formalism, however, these 
data in general belong to different analytic sheets 
of the two- variable Green's function (Fig. 12b I, 
due to the presence of an infinite set of branch 
cuts. The analytic sheets are defined on so-called 
wedges in the complex variable space C". Depend- 
ing on the considered wedge, an analytic continu- 
ation within the sheet to real frequencies and volt- 



3 



data from 

QMC 
simulations 



section , III, 



function theory of several complex variables 



set of 
wedges 



comprised 



give rise to 



l^ranch cuts of 
G in 



separated 
by (Fig. 







discrete data 


set G(ii/^,y,, iui-n) 




defines 




(appendix 




function 







of 



wedges 



Vladimirov's 
formula (43] 




theory of 
integral 
representations 



can be 
thought 
of as 



domains of 
holomorphy 



represents\ 



Bergman-Shilov 
boundary 



domains 


[Cauchy] 


(topological) 




represents 


boundary 



conventional function theory 



sectiQn ,IV, 



,viek|^ 



anaktic continiratiori 
Matsubara voltage theory 



real-time 
structure 
(set of edges 



continuity 
assu mptio n 
Eq. (gsI 



[MaxEnt 



approx. 


inferred 




probability 




distribution 



spectral function 
A(lo) 



transport 
properties 



section Jl, 



FIG. 1: Contents of the paper as a flowchart. Full arrows represent numerical procedures within 
our approach. Empty arrows represent the formal steps necessary for a derivation of the MaxEnt 
kernel used for numerical results. Shaded areas are discussed in the corresponding sections. The rather 
unconventional two-dimensional function theory with respect to the simultaneously variables {zip,Zu)) 
is explained by analogies (=) with the conventional theory in the light-shaded area. It may be skipped 
at first reading. Concepts such as the biholomorphic maps are used to obtain explicit equations. Dash- 
dotted lines denoting either Vladimirov's formula or Cauchy's integral equation put an emphasize on the 
fact that the content of a domain (of holomorphy) is parametrized by its (Bergman-Shilov) boundary. 
In the numerical implementation it is thus an inverse problem to reconstruct function values on the 
(Bergman-Shilov) boundary which requires a MaxEnt approach. 



ages may not have a direct physical interpretation. 
However, QMC data from these sheets should be 
used to reconstruct the physical real-time limit. 
For this sake, the mathematical structure of the 
Green's function's sheets is systematically inves- 
tigated in section [ITTj Since it may require a lot of 
effort to study the details of the latter, the reader 
is recommended to first study the results section 
and possibly skip section 1111] at first reading. 



fines the inverse problem ( 72 ) . It is shown in sec- 



The results section |TT] discusses numerical re- 
sults for the dot electron spectral function out of 
equilibrium and consequent transport characteris- 
tics. It relies on a MaxEnt procedure which infers 
a probability distribution based on a linear rela- 



tion ( 72 1 which is derived in section IV but has 



to be inverted. Since as in the conventional Wick 
rotation of imaginary-time data,'^^ the inversion 
is an ill-posed problem, the MaxEnt provides a 
most probable solution by means of Bayes' the- 
orem. Central to the function-theoretical deriva- 



tion of the relation ( 72 1 is a continuity assumption 



to the real-time structure of the theory, Eq. ( 65 ) 



It gives rise to the kernel operator ( 73 1 which de- 



tion |IV| that the continuity assumption improves 
the MaxEnt algorithm for the determination of 
spectral functions dramatically, as compared to 
the earlier approach introduced in Ref. I10[ such 
that nontrivial nonequilibrium spectra could be 
obtained. 

Let us now briefly discuss the more mathemati- 
cally involved section|TTl]by means of the flowchart 
in figure [T] Since the analytical derivation makes 
use of two-dimensional complex analysis, it shall 
provide an introduction to that field, comparing 
its fundamental notions to those of conventional 
function theory. In particular, the theory of in- 
tegral representations of functions on domains of 
holomorphy is discussed. Such functions may of- 
ten be parametrized by their values on the so- 
called Bergman-Shilov boundary. In an analogous 
way, conventional complex analysis parametrizes 
functions on domains by their boundary values 
using, e.g., Cauchy's integral equation. Also the 
concept of the conformal map has the analogon 
of a biholomorphic map which is widely used for 
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FIG. 2: Effective-equilibrium data as obtained 
from CT-QMC simulationpl for U = ST, e$ = 
O.ir, l3 = 5r~^. The integer number m speci- 
fies the respective index of the Matsubara voltage 
iPfn — Airm/lS. The discontinuity at a;„ — ipm/^ is 
the principal branch cut which will in particular 



be discussed in section III D 



will be discussed in detail in section IV and ap- 
pendix|Ej we here assume a certain property of the 
Green's function, namely that its real part at the 
meeting point of the cones is independent of the 
cone it was approached from. This allows to map 
data in different cones by means of linear transfor- 
mations into the actual data space and hence im- 
prove the accuracy of the MaxEnt tremendously. 
We must emphasize that the validity of this cru- 
cial property cannot proven rigorously; however, 
the results obtained can be taken as evidence that 
its violation does not influence the physical struc- 
tures too much. Furthermore, in section [V] we 
provide a route to improve on this approximation 
systematically, at the expense of a more complex 
algorithm. 

We concentrate on the evolution of the spec- 
tral functions as function of Coulomb parameter 
U and bias voltage $. Using the relation ^ we 
also calculate I{V) characteristics and compare 
them to results obtained with other techniques. 



formal derivations in the present work. Section 
in] also systematically points out that wedges are 
the domains of holomorphy for the two-variable 
Green's function and provides constraints to the 
Green's function which give rise to Vladimirov's 



integral representation (43) which is central to 
the constructed MaxEnt a gorithm. The repre- 
sentation links function values on the real-time 
boundary (the edge) of a considered wedge to data 
within the wedge. The dash-dotted lines in figure 
[l] represent such linear relations which are practi- 
cally used in the reverse direction and thus bring 
along an inverse problem. 



II. RESULTS FOR SPECTRAL 
FUNCTIONS AND TRANSPORT 



A. Weak-Coupling Regime 

The resulting nonequilibrium spectral function 
for U ^ 2T, P ^ 5r-\e$ = F, obtained by 
evaluating A{uj) = where A is some 

two-dimensional MaxEnt-inferred quantity, is dis- 
played in figure [Sj A good agreement with the 
zero-temperature second-order perturbation the- 
ory provided by Ref . [Sni is observed. Presumably 
due to the finite temperature, the quasi-particle 
weight is slightly smaller than the perturbative 
one. Also included is the default model sup- 
plied to the MaxEnt, a Lorentzian whose width 
is determined by the Bayesian procedure outlined 
in appendix [Ej For the weak coupling case, the 
width is essentially the one for the true spectrum, 
as expected. 



In the following, we present results obtained 
from CT-QMC data using the algorithm described 
in Ref. |10i Figure [2] provides an example for 
the raw simulation output. It has to be analyt- 
ically continued with respect to both, Matsub- 
ara voltage and frequency. This was accom- 
plished using MaxEnt applied to the functional re- 
lation between Matsubara-domain data and spec- 
tral function developed in section IV As already 



discussed in Ref. HOl the Green's function is an- 
alytical in certain cones in the four-dimensional 
variable space. Previously, we used only the cone 
closest to iifm = for providing information to 
the MaxEnt. This turned out to be not sufficient 



B. Intermediate-Coupling Regime 

Figure |4] shows spectral functions computed at 
U = 4r and inverse temperature /3 = 5r~^ for 
different bias voltages $. An excellent agreement 
with second-order perturbation theory is obtained 
for the cases e$ = 0.25r, e$ = O.SF, e$ F. 
Although one should not expect it to be very dif- 
ferent, the quasiparticle resonance in the MaxEnt 
result for e<I> = 0.125F is significantly reduced as 
compared to e<I> = 0.25r and disagrees with the 
almost unchanged second-order perturbation the- 
ory. This is a systematic MaxEnt artifact which 
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FIG. 3: (color online) Spectral function of the 
dot electrons as inferred for the nonequilibrium 
weak-coupling case U = 2T, e<i> = F, /? = SF"-"^, 
compared to zero-temperature second-order per- 
turbation theory. 



prefers to reproduce the default model in case of 
missing Bayesian evidence. In our case, Bayesian 
evidence is indeed decreased for very small volt- 
ages, because the simulation data at iprn — can- 
not be taken into account as the Matsubara volt- 
age Green's function has a branch cut there (see 
Ref. |5|for details). 

The normalization of the MaxEnt spectra is 
reasonably close to one. In particular at the small 
voltages, and for larger frequencies some side- 
bands form which tend to increase the total spec- 
tral weight unphysically. However, the good de- 
scription of the low-energy physics seems to be un- 
affected by this type of artifact. Note that again 
the optimal width of the default model is the same 
as the one for the final spectrum. 

Inferred spectral functions for an even stronger 
interaction U — QT are displayed in Fig. [5] for in- 
termediate to large bias voltage. For the equilib- 
rium situation we already are in a regime with a 
distinct three-peak structure with an Abrikosov- 
Suhl resonance (ASR) at w = characteristic for 
the Kondo regime. An estimate for the equi- 
librium Kondo scale gives Tk ~ O.IF. Again, 
the general low-energy behavior agrees very well 
with the results from perturbation theory - which 
should still be valid for this value oi U - and the 
weight at e$ = 0.125F is again underestimated. 
At e4> = 0.5F one now observes a distinctly larger 
broadening of the ASR as compared to perturba- 
tion theory, and at = F a clear double-peak 
structure is visible. This structure is compatible 
with a Kondo peak splitting. However, due to the 
approximations involved, we feel unable to decide 



at present whether this feature is actually a pre- 
diction of the Matsubara-voltage theory itself. It 
is interesting to note that here the default model 
is strongly renormalized for small bias, while at 
large bias the default model, apart from the dou- 
ble peak structure around cj = 0, again is already 
a reasonable estimate for the full spectrum. 

As before, the spectral weight in the now de- 
veloping Hubbard bands is strongly enhanced as 
compared to perturbation theory, pointing to- 
wards an overestimation of the integral weight by 
MaxEnt. On the other hand, the position is in 
good agreement with perturbation theory. 

Figure |6] shows a similar set of curves for in- 
teraction strength U = 8F. As compared to the 
lower values of U, a similar behaviour of the al- 
gorithm is observed. The perturbative prediction 
for the ASR is again essentially reproduced. The 
ASR is however again broadened as a function 
of the bias voltage and appears to split eventu- 
ally at e$ = F. Again, at the small bias voltage 
e$ = O.IF, lack of Bayesian evidence causes an 
underestimation of the ASR. The increased inter- 
action again broadens the spectra as compared to 
smaller values of U. However, the MaxEnt proce- 
dure does not clearly predict the correct Hubbard 
peak positions at zLU/2. A possible reason for this 
is the partially rather slow decay behaviour of the 
kernel function which was derived for the MaxEnt. 
It may result in decreased resolution in the high- 
frequency range, as compared to a conventional 
MaxEnt procedure for the Wick rotation. 

Spectra for interaction strength U = lOF are 
displayed in figure [7j Again, the solution is very 
similar to the perturbative prediction, and still 
a splitting of the ASR at larger bias voltages is 
observed. The MaxEnt resolution issue for the 
Hubbard bands is again observed. 



C. Approaching Lower Temperatures 

For data at lower temperature, namely /3 = 
lOF^^, the behaviour of MaxEnt solutions is sim- 
ilar to the one described above. Nevertheless, 
sharper structures, such as the Kondo peak, which 
emerge at lower temperatures, make the MaxEnt 
procedure more challenging. A common way to 
deal with this problem is the so-called "annealing 
procedure" .'2S»29j Here, a fine temperature grid is 
imposed in order to freeze out low-energy features 
step by step. The procedure starts with a feature- 
less default model at very high temperatures. At 
lower temperature, the MaxEnt result of the next 
higher temperature is used as default model, and 
so forth, until the target temperature is reached. 
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We found earlier'^ that it is of great use also 
within the two-dimensional analytic continuation 
problem. Also in the present extension of the ap- 
proach in Ref. [101 the MaxEnt yields more well- 
behaved solutions if a higher temperature is used 
as default model. This was investigated by a sim- 
ple single-step annealing procedure, using results 
from a. (3 — 5r^^ run. In fact, the occurance of 
unphysical normalization- violating sidebands may 
already be avoided in some cases for this rather 
rough temperature grid. Figure [8] shows two ex- 
amples in which the one-step annealing procedure 
was able to improve the results significantly. 



D. Transport Properties 

Using the Meir-Wingreen equation ([2]), we are 
able to compute transport properties based on 
spectral functions resulting from the MaxEnt an- 
alytic continuation procedure. Figure [9] compares 
results obtained at Coulomb interaction strengths 
U — AT and U = GF to real-time quantum Monte- 
Carlo data from Ref. 

The low-temperature data were obtained with- 
out the employment of an expensive annealing 
procedure: the most distinct feature of the low- 
temperature data is an increase in current at low- 
voltages for j3 — 10F~^ and U = QT. It is obtained 
no matter whether the single-step annealing pro- 
cedure described above is employed or not. 

As compared to the real-time QMC calcula- 
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tions, we obtain a good agreement at interaction 
strength U = AT. Since the only significant devia- 
tion is at e4> = 0.25, the discrepancy may still be 
due to statistical errors within the Monte-Carlo 
procedures. At higher interaction U = 6T, the 
current predicted by our method appears to be 
systematically higher at voltages e<i> > 0.2r, for 
both temperatures. 

In the following, we will discuss the low-voltage 
region of both cases, U — AT and U = BF. As 
shown in the previous section, in the case f5T ~ 5, 
the quasi-particle weight is underestimated in the 
voltage range e$ < 0.125r for the reason that 
i'/'m = data cannot be taken into account for the 
analytic continuation. This is compatible with the 
fact that the current is underestimated as com- 
pared to the (3T = 10 data, as well as the real-time 



QMC data. 

Figure [To] shows current- voltage curves for dif- 
ferent values of U at /3T = 10. In any case, 
an S-shaped current-voltage characteristics is ob- 
tained, the first increase of which is due to the 
ASR, and the second of which is due to the Hub- 
bard bands. The reduced weight of the ASR at 
increased interaction brings about an earlier de- 
parture of the U = ST curve as compared to the 
curves at U = 4r and U — QT. Same is true 
for the U — 6T curve as compared to the case 

U = 4r. 

In the next section of the paper, mathematically 
involved details of the structure of the Green's 
function will be analyzed systematically. It comes 
along with an introduction of the relevant con- 
cepts of the theory of several complex variables. 
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Those are in particular needed for the derivation 
of the MaxEnt inverse problem, equation (72), 



which leads to the results discussed in the present 
section. The latter mathematical derivation is 
provided in section |IV[ It may be convenient to 
skip section [ni] on first reading and directly go to 
the MaxEnt derivations and discussions of section 

HYi 



III. SEVERAL COMPLEX VARIABLES 
AND THE GREEN'S FUNCTION 

We want to provide a detailed physicist's intro- 
duction to the basic mathematical structures of 
functions with several complex variables, essential 
for a full understanding of the Green's function 



within the Matsubara-voltage formalism. In par- 
ticular, we will attempt to answer the question by 
what means intuition from conventional function 
theory is appropriate or misguiding in the context 
of dynamic expectation values. 



As next step we will describe the ana- 
lytic structure of dynamic correlation functions 
G( ) appearing as a fundamental object in 
the Matsubara-voltage formalism (see I for de- 
tails) by means of the function theory for multiple 
complex variables. Uniqueness of the mathemat- 
ical procedure involved in the analytical continu- 
ations is proven and the connection to maximum 
entropy approaches for the inference of spectral 
functions is made. 
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FIG. 7: (color online) Nonequilibrium spectral functions for U = lOF at inverse temperature f3 = bT ^, 
as compared to zero-temperature second-order perturbation theory. Line legends are the same as in 

Fig. [3} 
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(a) bias voltage e<I> = O.SF (b) bias voltage = l.Or 

FIG. 8: (color online) Lower-temperature spectral functions as inferred for U = 4F, /3 = lOF^^ with or 
without an annealing step. The default model for the single-step procedure is taken from temperature 

/3 = 5F-1. 



Holomorphy of Functions of Several 
Complex Variables 



Rieniann differential equation 



In order to discuss the notion of holomorphy 
in the context of functions with several complex 
variables we will partially follow the book by 
Vladimirov.l^ 



A function G of one complex variable is holo- 
morphic at a point zq, if and only if the Cauchy- 



dG 



= 



(3) 



is satisfied. The notion for holomorphy of func- 
tions of several complex variables is a natural ex- 
tension of this definition: A function /(z), with 
z € C*, is holomorphic with respect to z at the 
point z^*^' if and only if it is holomorphic with 
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FIG. 9: (color online) Transport characteristics as 
compared to real-time QMC data from Ref. US) 



T 




Voltage e<I>/r 

FIG. 10: (color online) Current as a function of 
the voltage for different interaction strengths. 



respect to each individual variable, 



Note that in the following, we will always denote 
vectors in <C^{d > 1) by an underlined symbol 
such as z. Hartogs ' fundamental theorem asserts 
that definition Q is also equivalent to the Weier- 
strafi definition of holomorphy for several vari- 
ables. The latter calls a function holomorphic 
at Zq if and only if there exists an open neigh- 
bourhood A/ of 0^°^ , such that for all z € M the 
function / may be written as an absolutely con- 
vergent power series f{z) = X]|^|=o ~ -^°'')"- 
a denotes the multi-index for the monomial z" := 

n^=o -^n"' ^^'^ 1*^1 ■= J2t=Q An analytic com- 
plex function of several variables is holomorphic. 
The notion of holomorphy implied by Eqs. Q 



is thus as natural and intuitive as in the one- 
dimensional case. 



B. Domains of Holomorphy and 
Biholomorphic Transformations 

The major qualitative difference between 
single- and multi- variable complex analysis is con- 
tained in the notion of a domain, and the geo- 
metric equivalence among holomorphic functions 
arising from classes of domains. This has far- 
reaching consequences to the theory itself, such 
as the construction of integral representations. In 
the context of our formalism we aim at integral 
representations. We will thus first comment on 
the basic structures which integral representations 
operate on. Furthermore, we will point out the 
most prominent differences to conventional func- 
tion theory. The notion of a domain in single- 
variable complex analysis is replaced by the no- 
tion of a domain of holomorphy in multi-variable 
complex analysis - the notion of a conformal map 
is replaced by the notion of a biholomorphic map. 
We will address the most prominent differences by 
first reminding the reader of basic structural prop- 
erties of one-dimensional complex analysis and 
then introducing the corresponding terminologies 
of the more-dimensional theory. 



1. One-dimensional function theory 

In one-dimensional complex analysis, domains 
are defined as open connected subsets of C. For 
the time being, we will restrict the discussion to 
simply connected open sets, i.e. open connected 
sets with no holes. 

Conformal maps between domains U, V, 
namely functions 

m : U ^ V, z^ m{z) (5) 

which are holomorphic and invertible (one-to- 
one), provide links between certain classes of do- 
mains. The Riemann mapping theorem states 
that conformal maps between any simply con- 
nected ?7 7^ C,0 and the unit disk exists. I.e. all 
simply connected domains 7^ C, are conformally 
equivalent: their structures of holomorphic func- 
tions map one-to-one to each other and are con- 
formally diffeomorphic. Generalizing the concept 
of holomorphy to Riemannian surfaces, conformal 
maps exist only for surfaces of the same topologi- 
cal genus. The uniformization theorem finds that 
for simply connected Riemann surfaces (topolog- 
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ical genus 0) up to conformal equivalence three 
classes exist: 

• the unit disk Di = {z e C : \z\ < 1}, 

• the complex plane C, 

• the Riemann sphere C U {oo}. 

These three surfaces form the so-called moduli 
space of genus 0, defined as the space of confor- 
mally inequivalent Riemann surfaces of genus 0. 

In general, the size of the moduli space of a 
Riemann surface of genus g grows as a function of 
g. Each modulus represents an equivalence class 
of holomorphic functions. 

Conformal equivalence plays an important role 
in physical applications such as two-dimensional 
potential flows around airfoils or conformal quan- 
tum field theory. 

One natural consequence of the conformal 
equivalence of all non-empty simply connected do- 
mains [/ C is that there always exists a func- 
tion which is not analytically continuable beyond 
the domain: the function 

CO 

/o(z) = ^z"' (6) 

a=0 

is holomorphic on the unit disk but may not be 
analytically continued to larger domains.'^ Using 
a conformal map m from the unit disk to U, which 
exists due to the Riemann mapping theorem, one 
finds the function /o om~^ which cannot be ana- 
lytically continued beyond U . One also calls the 
unit disk the domain of holomorphy of fo, i.e. the 
largest domain for which /o is holomorphic. U 
is the domain of holomorphy of fo o m^^. For a 
given domain G, if there exists any function / for 
which G is the largest possible domain in which / 
is holomorphic, the domain is called a domain of 
holomorphy. 

In general, any domain is also a domain of holo- 
morphy in conventional complex analysis. Due 
to the simple structures arising from these far- 
reaching equivalences, the conventional function 
theory of one complex variable is widely consid- 
ered to be a finalized field of mathematical re- 
search. 



2. Multi- dimensional function theory 

The two final statements of the last section are 
completely incorrect for several complex variables 
C^, d>2. 

In several complex variables. 



• a domain is not necessarily a domain of 
holomorphy; 

• domains of holomorphy are usually not bi- 
holomorphically equivalent. 

An example of a domain which is no domain of 
holomorphy is the hollow sphere (see paragraph b 
below), because all analytic functions in the hol- 
low sphere can be analytically continued to the 
sphere. As compared to one-dimensional func- 
tion theory, the emerging structures are thus very 
rich. Depending on the domain geometrjy, very 
different sheaves of holomorphic functions^^ will 
arise. 

a. Biholomorphic maps. The tool of a bi- 
holomorphic map, as mentioned in the second 
point, is the generalization of a conformal map to 
several complex variables. If a holomorphic map- 
ping m : f/ — )• [/', with U, U' C is invertible, it 
is called a biholomorphic map. Two domains C/, 
U' are biholomorphically equivalent, if and only if 
such a map exists. They will have an equivalent 
sheaf of holomorphic functions. Biholomorphic 
maps do not necessarily preserve angles. There- 
fore, they are usually not conformal. Neverthe- 
less, with respect to the holomorphic structure, 
they are the natural generalization of conformal 
maps on C, because they are complex diffeomor- 
phisms: the inverse of a biholomorphic map is also 
holomorphic ^ 

b. Domains of holomorphy and holomorphic 
envelopes. A striking example of a domain C 
which is not a domain of holomorphy is given by 
the hollow sphere M := {z = (^1,^2)^ G : 
i < ^l^ip + |z2p < !}■ In sharp contrast to the 
single-variable case, one can show that any holo- 
morphic function / : M — > C may be analytically 
continued to the unit sphere. The unit sphere 
is, in fact, a domain of holomorphy and is thus 
named the holomorphic envelope of M.^i^lxhis ex- 
tends to the famous result by Friedrich Hartogs 
that isolated singularities are always removable 
for analytic functions — ^ C, d > 2 {Hartogs' 
lemma). While isolated singularities play an es- 
sential role in the residue calculus in the d = 1 
case, the d> 2 case is, due to Hartogs' result, of 
an entirely different nature. As we will see in sec- 
tion [IIIC] the theory of integral representations of 
complex functions for d > 2 has consequently a 
very different character as compared to the d = I 
case. 

Since ordinary domains such as the hollow 
sphere are rather friendly as far as analytic con- 
tinuation is concerned, mathematicians restrict to 
the systematic study of corresponding envelopes 
of holomorphy, i.e. the domains of holomorphy. 
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One can show that a domain is a domain of holo- 
morphy if and only if it is a so-called pseudocon- 
vex domain. Pseudo convexity is a certain gener- 
alization of convexity from to C^^l^^l For this 
reason, the fundamental domains of holomorphy 
in our application, the wedges, are indeed convex 
(cf. section III D ). 

c. Biholomorphic equivalence. In order to 
provide a classification of domains of holomor- 
phy in C^, d > 2, the fruitful strategy from 
d — 1, namely using biholomorphic equivalences, 
is adopted. Non-empty simply connected domains 
of holomorphy C are usually not biholomor- 
phically equivalent. There are many different 
types of holomorphic structures, depending on the 
domain geometry. 

A prominent example of biholomorphically in- 
equivalent domains of holomorphy is given by the 
unit ball \z\ < 1 and the bicylinder Di x Di in C^. 
The unit ball in higher dimensions is with regard 
to the holomorphic structure in no way related to 
the unit disks and carries in fact a rather singular 
holomorphic structure. 



C. Integral Representations 

As one might have already guessed from Har- 
togs' lemma and the domain dependence of the 
mathematical structures, finding an analogy to 
Cauchy-like integral representations, yielding, for 
example, the spectral representation of our two- 
variable Green's function, turned out to be a cum- 
bersome task. 

At present, one knows several, more or less 
general, possibilities to construct such represen- 
tations, even with different dimensions of the 
integration manifolds. The Bochner-Martinelli 
representatiorP^ is, for example, probably the 
most general integral representation, but the in- 
tegration manifold is 2(i — 1 dimensional. 

For our practical purposes, a minimal integra- 
tion space dimension is, of course, most desir- 
able to reduce the number of fitting variables 
when reconstructing the function by a Maximum- 
Entropy-like Bayesian inference technique. 



1. The Shilov Boundary 

In d = 1 complex analysis, Cauchy's integral 
representation reconstructs all values of a simply 
connected open domain of finite radius using the 
limit values on the topological boundary, as long 
as the boundary values are continuous. I.e. the 




(a) Cauchy {d 



(b) Bergman-Weil (d -- 



FIG. 11: Comparison of Cauchy and Bergman- 
Weil integral representation theories. Integrations 
run over the (a) full topological boundary dD and 
(b) the Bcrgman-Shilov boundary S C dD of an 
analytic polyhedron, respectively. 



structure on the topological boundary determines 
the structure of the interior. 

For d > 1, the topological boundary can in 
this sense be reduced to an often much smaller 
set, the so-called distinguished or Bergman- Shilov 
boundary]^ln the 1930s, Stefan Bergman discov- 
ered the distinguished boundary in the context of 
for bounded domains with piecewise smooth 
boundaries. He found that under certain regular- 
ity conditions, an integral representation with re- 
spect to an only two-dimensional manifold S of all 
intersections of the smooth boundaries was possi- 
ble, using the so-called Bergman kernel function. 
A precise geometric and analytic description of 
these results is provided in great detail in chapter 
XI of Ref. HH A pictorial comparison of the result- 
ing Bergman-Weil representation to the Cauchy 
formula is provided in figure |11| These concepts 
where independently discovered by Shilov in the 
1940s in the rather different context of commuta- 
tive Banach algebras.'^ 

a. Generalization. In the modern terminol- 
ogy, the Shilov boundary, as a generalization of 
the Bergman-Shilov boundary, may be defined 
for any compact space with respect to an alge- 
bra of continuous complex-valued functions on 
the space.l^ If we for example find physical con- 
straints imposing certain conditions to the set of 
holomorphic functions, the Shilov boundary with 
respect to these functions may be reduced to a 
smaller set. If a sufficiently elaborate kind of in- 
tegral representation is used, this may enable us to 
again reduce the number of linear fit parameters 
significantly in the Bayesian inference problem. 

b. Examples. For the domains with piece- 
wise smooth boundaries which Bergman investi- 
gated, S is given by the unification of all pos- 
sible intersections between the smooth boundary 
hypersurfaces, as long as certain regularity condi- 
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tions hold. We return to our two examples from variables itpm, iw„ readi^' 



part III B 2 c the bicylinder and the unit sphere 
and comment on their Bergman- Shilov bound- 
aries. 

The bicylinder, DiX Di, is one of the most eas- 
ily accessible domains of holomorphy, because it 
simply factorizes into two Di-disks in C. A mini- 
mal integral representation of a holomorphic func- 
tion which is continuous on the closure Di x Di, 
is simply given by the product of two conven- 
tional Cauchy representations (see Theorem 2.2.1 
in Ref.lH]), 



1 



(27ri)2 



/(C) 



Six Si 



Uki^k-Ck)' 



(7) 



Therefore, the Bergman-Shilov boundary of Di x 
-Di is given by the only two-dimensional toroidal 
subset S'l X 5*1 = dDi x dDi of the three- 
dimensional topological boundary d{Di x Di) — 
dDi X Di U Di X dDi . Similarly, integral repre- 
sentations of domains which are direct products 
of domains can be constructed from the con- 
ventional Cauchy integral formula. 

It was already mentioned in part |III B 2 c| that 
the unit sphere in is not biholomorphically 
equivalent to the bicylinder. In contrast to the 
two connected smooth boundary hypersurfaces of 
the bicylinder, the boundary hypersurface of the 
unit sphere is not even a smooth hypersurface in 
Bergman's notion. One can show that here the 
Bergman-Shilov boundary is, in fact, identical to 
the topological boundary. Thus, any integral rep- 
resentation for holomorphic functions in domains 
such as the unit sphere must invoke at least three 
real integrals. It is an example of the strong dis- 
tinguishments which have to be made between 
certain classes of domains of holomorphy. 



D. Holomorphic Structure of the Green's 
Functions 

As a next step, we will systematically analyze 
the mathematical structure of the Green's func- 
tion arising in the Matsubara-voltage approach. 
The bare Green's function with respect to the two 



E 



(8) 



^■^^ iw„ - a{iip,n - $)/2 - £d + irsgn„ 



Here, sgn„,„ := sgn(a;„ - aip,n/2). Performing 
the analytic continuations ia;„ z^^ iip„i — >■ z^, 
the sign function in the denominator results in an 
ambiguity, as far as the definition of domains, for 
which Go is holomorphic, is concerned. 

Choosing a branch cut structure which corre- 
sponds to the continuation 

sgn(a;„ ± <^to/2) sgn(Im ± Im z^/2) (9) 

appears to be most sensible from both, a mathe- 
matic and a physical point of view: 

From the former perspective, in contrast to 
other choices the resulting domains are also do- 
mains of holomorphy and are thus "maximal" 
with respect to the holomorphic structure. The 
four domains of holomorphy are given by sep- 
arated into wedges by the two branch-cut hyper- 



planes Imz^ ± ImZ;^/2 = (see Fig. 12a I. From 



the latter, it is just the imaginary part of the 



linear combinations of 



and 



appearing m 



the denominators of perturbative expressions in 
U, which yields the crucial sign-switching delta 
functions generating the non-analytic terms sep- 
arating the domains for which G is holomorphic. 
See, e.g., appendix B of paper I. 

Consequently, for a finite interaction U 0, the 
appearance of a branch cut must be expected for 
any new kind of linear combination of z^^ and z^p 
in the denominators of the integrands in the per- 
turbation expansion. For example, second-order 
perturbation theory from Ref. Ill indicated that 
particle-hole bubbles create higher-order branch 
cuts Imz^ — ^Imz^ = 0, with 7 € Z odd, due to 
the structure of convolutions which are involved. 
In a later publicatiorP^ it was pointed out that 
vertex corrections seem to introduce yet another 
branch cut for Imz^p = 0. The physical retarded 
Green's function is then given by 



limx .0+ + limA .0- 
G,et(^)= lim " . ^^^G($ + iJ^,z<,), (10) 

z„-!.w+iO+ 2 



i.e. one has to average over the two possible limits with respect to (5, 
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Note that the latter subtlety was not taken into 
account in the direct continuation using the cone 
Cg in Ref. IIOI and also not in the initial approach 
in Ref. E 



1. Resulting Mathematical Assumptions 

The following assumptions are being made for 
the mathematical structure of the Green's func- 
tion: 



1, 



By means of holomorphy, we obtain cone- 
like constraints for the combinations of 
as depicted in Fig 



12b 



imaginary parts 
More precisely, we require G to be a solu- 
tion of the Cauchy-Riemann equations Q 
for any z^^"^ except for those x^^-* for which 



lmz(0) 
Imz.fo) 



— Imzfi*' or 
2 ^ 

0, 



(11) 



with some 7 € 2Z+1. Those ^(o' define 
the branch cut hyperplanes and delimit the 
wedges for which G is holomorphic. 

2. We will require the interacting Green's func- 
tion G{z) to be bounded, i.e. 

sup |G(z)| < cx). (12) 

3. We will assume that the Green's function 
G{z) is uniquely defined by the discrete 
function values G(i(^m, i'^n) which are ob- 
tained from the effective-equilibrium com- 
putations. I.e. we require that the continua- 
tion to a multisheeted holomorphic function 

G(i(/3m,iw„) G{z 

) (13) 



IS unique. 

The second assumption is justified by the struc- 
ture of the convolution equations in perturbation 
theory and the boundary conditions that terms 
gj^^ Qiuinfs ^ evaluate to 1 and -1 before the 
analytic continuations are carried out. 

A proof of the third statement, which is of 
course crucial for the physical theory itself, will be 
provided in appendix [X] It is based on assump- 
tion 1 and 2 and the assumption 3' which sharpens 
the requirements on the z — >■ 00-asymptotics: 



3'. Given arbitrary x^^"^ 
we have 



e ]R2 \ {0} and C e 



lim 

f 00 



C-G(Cx(°)) 



< 00 



(14) 



In other words, G(z) is required to decay like 
a usual Green's function as a function of 
where z — (- x^^^ , if and only if x^^^ satisfies 
the regularity condition xi?) ^ ixlp"^ /2. 



2. Justification of assumption 3' 

The assumption 3' may be justified as follows. 

Consider the absolute value of the free Green's 
function (Isl), 
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\Go{z^,Zu,)\ < ^ 



a=±l 



a{z^ - $)/2 - Cd + iFsgnlm {z^ - az^/2)\ 



r 



It is obvious that it decays cx ^ when z^p 
and = C ■ for the nonsingular combinations 



of 



and 



It does not decay at all in the 



singular cases xi?^ — ±x)^' /2 



It is easy to check that interaction [/ > does 
not change this high-energy structure. Let us 
examine the second-order self-energy expression 
(Eq. (15) in Ref.[T]): 



E(2)(z) = [/2^ 



(ai - Q!2 + aa)- 



ei + 62 - £3 



r 



with fa, = /(e» - a,^/2). 

Due to the structure of the denominator, we see 
that on top of the singular directions of the bare 



Green's function, 



„(o) 



singular directions x^^ — ± 

Consequently, assumption 3' is incorrect for the 
second-order self-energy. Nevertheless, when in- 
serted into Dyson's equation. 



±x^^^ /2, we also have the 



-3^(0) 



G(2)(z) 



Go{z) 



l-GoU)S](2)(z)' 



(15) 



-3^(0) 



we see that for the directions xi 
limiting behaviour of Go is adopted, i.e. the be- 
haviour ( 14 ). 



Note that the uniqueness proof of the appendix 
also holds when directly applied to the self-energy, 
because the singular directions of the second-order 
perturbation theory, — ±^x^\ are not re- 
quired to be regular in the proof. This is because 



(0) 



the direction x, 
cut (assumption 1). 



i 2 



(0) 



also defines a branch 



E. Tubular Cone Domains ("Wedges") 

As we have seen, the structure of Go combined 
with the structure of convolutions in the pertur- 
bation theory with respect to U, indicates that 
the numerous branch cut hyperplanes divide 
into several, in fact infinitely many, wedges of the 
form = -f iG. G is by definition a convex 
cone with its vertex at zero. See also the picto- 
rial discussion in Ref. [TUl Due to the convexity 
of G, T*^ is pseudoconvex and thus a domain of 



holomorphy.^'^ In the mathematical classification 
scheme, domains like these are called tubular cone 
domains. 



1. Geometry of the Cones 

We briefiy introduce certain notions of the de- 
scription of the analytic geometry of cones in M.'^. 
This is necessary to thoroughly follow the math- 
ematical formulae which are involved in the de- 
scription of the analytic structure of T'^' . The 
cone G with vertex at zero is formally defined by 
the scaling property x £ C ^ VX > : Xx £ C . 

Its dual cone C* is defined via the standard 
scalar product 



(€'3^)-!]^^^'=; (16) 



fc=l 



by 



G* {e e M-^IVa; e G : (e,x) > 0}. (17) 



G* represents the space of positive semi-definite 
linear functionals on G when the functional form 
(16) is considered. The dual cone is important, 



because the construction of kernel functions often 
involves Fourier transforms. 

A convex cone is a cone for which the straight 
line between any pair of points within the cone is 
also contained by the cone. We will also use the 



analytic continuation of the scalar product ( 16 ) 
with respect to x. We continue x z holomor- 



16 



phically in ( 16 1 



M 



Z 



{i,z):=Y^^uZk; ie^^.^zeC". (18) 



k=l 



2. Analytic Structure and Biholomorphic 
Equivalence Classes 

Tubular cones and domains are well-known ob- 
jects in the theory of several complex variables, 
because they naturally arise in certain fields of 
mathematics. As a consequence, many efforts 
were put in for a detailed understanding of their 
structure. In the prominent physical example, 
axiomatic quantum field theory, the cones repre- 
sent forward and backward light cones, in four- 
dimensional spacetime, d ~ 4. A celebrated result 
was Bogolyubov's edge-of-the-wedge theorem.!^ 
Mathematical examples include Fourier analysis, 
functional analysis, the theory of hyperfunctions, 
and the theory of partial differential equations. 

A key component of the wedge is given by its 
edge, namely the real subspace associated with 
the vertex of the cone. Because in our case the 
vertex is located at zero, the edge of the wedge 
T*^ can be formally identified with an oriented 
copy of the real subspace. 



Edge-pc := 



iO^ 



(19) 



where O'-' is an infinitesimal vector within the cone 
C. Although there exist, depending on the direc- 
tion approaching the origin within C, several in- 
finitesimals 0*-^, the Edge is well-defined through 



(19), because all infinitesimals in C are obviously 



equivalent with respect to holomorphic continua- 
tion in T*^. T'^ may be regarded as a generaliza- 
tion of the upper half plane. 

Let C, C" be arbitrary convex cones in M'^, 
C ^ C. and T^' are in general not biholo- 
morphically equivalent. This means that the sets 
of holomorphic functions living on them are struc- 
tured differently. Wedges within the space are 
fortunately an exception to this rule: and T*-^ 
are biholomorphically equivalent for any combina- 
tion C',C'. See also the introductory notes in the 
corresponding part of the second volume of ref- 
erence TS! In C^, biholomorphisms between 
and T*^ niay be constructed easily using com- 
plexified rotations and dilations of E^: Consider 
that the real (non-singular) matrix M : — > 

]R^ ^ ^ ^ maps C to C", i.e. C" = MC. 

The corresponding biholomorphism between 
and is obtained from the complexified map 



z. It is easy to see 



that T^' = MT^, that M is holomorphic, and 
invertible. 

If the complexified linear map is a rotation, we 
will also call it a biholomorphic rotation when we 
want to emphasize the biholomorphic character of 
the mapping. 

Note that the helpful notion of the Bergman- 
Shilov boundary is not directly applicable to , 
because T*^ is unbounded. However, as will be dis- 
cussed next, sequences of bounded domains ap- 
proaching T*^ from its interior may be used to 
understand the holomorphic structure on . 



3. Bergman-Weil Representations 

A sequence of bounded domains 2?„ C T'^, 
lim„_j.oo 1^71 = , with piecewise smooth bound- 
aries like in Fig. |llb| may be easily constructed, 
such that the edge of the wedge contains a part 
of the Bergman-Shilov boundary 5„ of I?„ and 
such that the other subsets of 5„ disappear to 
oo as n — )• 00. In our case, d = 2, an ex- 
plicit construction of such a sequence may be ob- 
tained by intersecting T*^ with a growing bicylin- 
der Bn := i?„ • (£>! x Di), X>„ S„ n , with 
radius i?„ oc n. The procedure is sketched in 



Fig. 13 One finds that the sequence of Bergman- 
Shilov boundaries improperly converges to 



Sao Edgerpc U {oo}, 



(20) 



where "oo" shall informally denote the point or 
merely a list of points which emerge when 
is, depending on the holomorphic structure, com- 
pactified suitably. At first glance, each direction 
for approaching oo might yield a different point 
in {oo}. The points "oo" carry the additional in- 
formation which is necessary to turn the structure 
on EdgCjip into a unique description of the holo- 
morphic structure on . 

The Bergman representation for an analytic 
polyhedron in may be obtained by the follow- 
ing, rather technical, procedure, whose details are 
not particularly relevant but enable us to investi- 
gate rather explicitly the structure of a suitably 
large class of holomorphic functions on a wedge. 
An excellent pedestrian's introduction to it is pro- 
vided by Bergman's original monograph, Ref. 1141 
The book also provides a comprehensive introduc- 
tion to the Bergman-Shilov boundary and biholo- 
morphic maps, using the example of analytic poly- 
hedra. Equations of the form 



C''-) =p)(z,A,); 



(21) 
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FIG. 13: (color online) An artist's impression of 
asymptotically filling the wedge T'^ with the se- 
quence I?„. 2?„ is a domain with piecewise smooth 
boundaries created by intersecting with a grow- 
ing bicylinder i3„, creating the analytic polyhe- 
dra 2?„ . The corresponding sequence of Bergman- 
Shilov boundaries improperly converges to the set 
Edge U {00} . Therefore, the sheaf of holomorphic 
functions on T*^ is solely characterized by its edge 
values and its asymptotic behaviour 00. 



with Afc e Afc C M, fc = 1,. . . where ^^^\z,\k) 
are A-parametrized families of analytic functions 
of z, shall define the analytic polyhedron. Each 
equation yields a surface in z-space for a given 
A and a hypersurface hu in z-space as A is var- 
ied continuously. The mutual intersections Ski ~ 
hkf^hi yield the Bergman-Shilov boundary surface 
^ — Ufc I ^ki- Ski is then parametrized by a func- 
tion z — g^''^\\k,h)- A holomorphic function / 
on the analytic polyhedron may then be written 
with respect to the Bergman-Shilov boundary of 
the latter using the Bergman kernel function 



Bki — 



d{Xk, A( 



02 



fljzi, Z2, \l)fk{zi, \k) 

, (fl(-)-zO(Br^-.2) 
fkizi, Z2, Xk)fiizi, Qi''''\ Xl) 



- (22) 



2l)(0 



(feO 



Z2) 



The integral representation with respect to f\g 
then reads 



m 



1 

"8^ 



E 

k=^l 



(kl) (kl) 
02 



kl 



/A.a, f«(z,Afe)f(0(z,AO 
and can be applied directly to our 2?„ domains 



(23) 



An explicit test on whether a Bergman inte- 
gral representation for Go on the = S„ n 
domains is feasible yields that the subsets of 5„ 
which go to 00 may not be neglected for G{z) 
functions. This is because Go{z) has a nonzero 
limit as z 00 if one goes along the cross-shaped 
submanifold Hezi^ = ±Ilez^/2 and keeps Imz 
constant. Due to the independence of Imz this 
problem occurs for each of the wedges. Hence, 
the Bergman kernel function (22 1 is only of lim- 
ited use for us. We will thus not go into further 
details of this rather clumsy computation here. 

The formal use of the sequence I?„ enables 
us to see very explicitly that the Edge of T*^ is 
with respect to the representation ( 22 ) the only 



carrier of structural information which involves 
finite values of z. The rest of the information, 
then uniquely defining the holomorphic structure 
on T^ , is encoded in the several possible classes 
of limiting behaviour as z approaches infinity. 



4. Cauchy-Bochner Integral Representation 

As a straightforward consequence of this, as- 
suming a certain limiting behaviour of the con- 
sidered set of functions on T*^, integral represen- 



tations with respect to the Edge 



iO may 



be derived. Even more generally, a constraint on 
the function class which also limits the Edge be- 
haviour, can be imposed in such a way that the 
Edge function yields a unique description. The 
several possible z — 00 behaviours are then, us- 
ing the information from the constraint, encoded 
in the Edge. This appears to be linked deeply to 
the extension of the notion of the Bergman-Shilov 
boundary to the notion of the Shilov boundary 
mentioned in the course of section IIII C II Con- 
sidering a subset of the sheaf of holomorphic func- 
tions on a certain (compactified) domain may 
cause the Shilov boundary to "shrink" . 

We are now going to discuss one of the ear- 
liest developments going beyond simply restrict- 
ing the considered function set such that a naive 
extrapolation of Bergman representations, like in 
Fig. [l3j holds. It was another extension of the 
Cauchy integral formula with respect to tubular 
cones by Salomon Bochner. He considers a func- 
tion / which is holomorphic on T*^ and satisfies 
the constraint 

11/(3; + iy) II <M,,;(C")e^l^l (24) 

which has to hold for any compact cone C" C C 
and for any e > 0, where M^j{C') is a suit- 
ably chosen real number. || • || is a norm which 
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integrates out the x variable, \\f{x + iy)|P := 
/d'^x|/(x + i7;)|2P 

Inequality ( |24[ ) constrains the limiting be- 
haviour z -> oo in a sufficiently strong way such 
that an integral representation with respect to the 
Edge may be constructed. 

Namely, the Cauchy-Bochner representation 
then allows a function / satisfying ( 24 1 to be writ- 
ten as 



/U) = 



1 



d'*a;'A:cU-a^')/(x' + iO^); (25) 



with the Edge values /(x' -f iO*^). Here, the so- 
called Cauchy kernei^'^ of the cone C, defined as 



(26) 



was introduced. It is straightforward to compute 
the Cauchy kernel for our wedges with this for- 
mula: see the Appendix of Ref. [TOl For our pur- 
poses, we will provide a general but easily applica- 
ble expression for further numerical and analytical 
computations in section [ill G| 

Unfortunately, as in the Bergman approach, a 
numerical test of Eq. ( 25 1 for / ~ Gq using an ar- 



bitrary wedge for which Go is holomorphic shows 
that the Cauchy-Bochner representation (25 1 is 



also incorrect for Gq. As a consequence, wc find 
that the Green's functions does not satisfy (24). 



This is compatible with the fact that the left hand 
side of ( 24 ) diverges in the case / = Gg, no matter 
which G' is considered. 

Nevertheless, as we will see, the Cauchy- 
Bochner kernel ( 26 1 will serve as a building 



block for the construction of an exact integral 
representation for a different class of holomorphic 
functions which in fact contains the holomorphic 
branches of our Green's function G{z) on the re- 
spective wedges. As such it is essential as a con- 
nection of real-time and imaginary-time structure 
of the Green's function. 



5. The tubular octant EI X EI and 
Biholomorphic Equivalence to the Bicyhnder 

Due to the biholomorphic equivalence of all 
in any of our wedges which arise for the Green's 
function may be mapped biholomorphically to the 
tubular octant H x H = i (K+ x K+), where 
H is the upper half plane of C. This domain may 
itself be mapped biholomorphically to the bicylin- 
der Di x Di via a piecewise Mobius transforma- 
tion of the coordinates. Hence, all wedges of the 
Green's function are biholomorphically equivalent 




mc.(Edge 



Si X Si 



FIG. 14: On the distinguished boundary surface 
of the bicylinder, the picture of Edge^-c under 
the component-wise Mobius transformation mc 
is delimited by the dash-dotted lines poo- A dis- 
continuity of G\j,c o Tn^^ occurs at the intersec- 
tion point (1, 1) of the two circles and prevents 
a Cauchy representation from being applicable. 
The violation of the Cauchy-Bochner condition 
( 24 ) appears to be related to the occurance of the 
discontinuity. 



to the bicylinder. Let us denote a corresponding 
biholomorphism by mc : — > I?i x £>i. We 
would like to comment on this due to the striking 
simplicity of the bicylinder and of domains which 
are direct products of domains with respect to 
the construction of integral representations. 

From the point of view of this construction, 
the tubular octant may be regarded as the sim- 
plest representant of the biholomorphic equiva- 
lence class of all wedges in C^. 

Due to the premises of the Cauchy integral for- 
mula, a usage of the representation ([7| for the 
biholomorphically transformed sheaf of holomor- 
phic functions is feasible in case the transformed 
Green's function G\j,c o tn^^ is continuous on its 
topological boundary 9(Di x Di). Note that un- 
der the biholomorphic transformations, 5*1 x Si is 
mapped to the edges of the T'^ wedges. The "oo" 
in maps to "oo" in H x H under biholomorphic 
transformation, and that again maps biholomor- 
phically, using the component-wise Mobius trans- 
formation 1^ to the points poo := {1} x Si U 
5*1 X {1}. Note that the points poo C d{Di x Di) 
are also part of the distinguished boundary torus 
X Si. 

The boundary behaviour of the transformed 
Green's function G\j,c °^c^ is not continuous at 
the intersection point (1, 1) of the two circles poo, 
due to the properties of Golj^c at oo leading to 
singular directions, as summarized in assumption 
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3' of section \IlTlJ\ 

An illustration may be found in Fig. 

Therefore, using the biholomorphic equivalence 
to the bicylinder is not immediately helpful for the 
construction of an integral representation of the 
Green's function G. Nevertheless, it is essential 
in the application of Vladimirov's approach which 
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will be subject of the next section. 



F. Vladimirov's Integral Formula 

Vladimirov provided a generalization of the so- 
called Herglotz-Nevanlinna representation for the 
upper half plane to tubular cone domains. His 
investigations were motivated by applications in 
the field o f linear passive systems in mathemati- 
cal physicj^ ° * ^'^ * ^° l Due to its generality, the ap- 
proach is applicable to the analytic wedges of the 
interacting Green's function G. 



1. Herglotz-Nevanlinna Representation (d=l) 

We will discuss the conventional Herglotz for- 
mulae. Herglotz' representation theorem consid- 
ers holomorphic functions in the open unit disk 
Di which have a positive real part, the so-called 
Caratheodory functions^-^*'^. By separating a 
phase factor out of the function one can also con- 
sider functions with positive or negative imagi- 
nary part, and so on. Since the open unit disk can 
be conformally mapped to the upper half plane 
H, using the Mobius transformation as men- 

' ^ Z + l ' 

tioned above, the representation can under certain 
circumstances be also used for H. 

By considering Caratheodory functions, Her- 
glotz' theorem only imposes assumptions on the 
positivity of the real (imaginary) part of the 
function. In contrast to Cauchy's integral for- 
mula, no assumptions about the behaviour of the 
Caratheodory functions on the boundary of the 
disk are made, such as the continuity. 

The theorem states that every Caratheodory 
function / can be represented by 



A comparison of Eqs. (27) and (28) shows that 



/(z) = i-Im/(0) + 



2" ei* + z 



Mt) , (27) 



where dcr is a nonnegative finite measurc^^^ 

Regarding a different set of functions, a for- 
mally very similar representation is the so-called 
Poisson formula, which is the analog of the 
Cauchy formula to the real analysis of harmonic 
functions (solutions of Laplace's equation), and 
can in fact be derived from it. It provides an inte- 
gral kernel for the solution to the Dirichlet prob- 
lem for the Laplace equation on the unit disk in 
M^. For a continuous function f : dDi — > M it al- 
lows to construct a harmonic function u : Di — > M 
as follows (see pp. 169ff. in Ref . ITS]) : 



, , 1 f^'' e'* -I- z 



ait 



z 



f(e")dt 



(28) 



the measure dcr of the Herglotz formula is in fact 
defined by the (possibly singular) boundary limit 
of the holomorphic function. 

As a natural extension of the Poisson formula, 
the Schwarz integral formula reconstructs a holo- 
morphic function / on the closed unit disk from 
the real part of its boundary values, up to a con- 
stant imaginary offset. It reads (p. 171 in Ref. fT5|) 



/(z)-iTm/(0)- 



1 

'2tt 



e'* + z 



/(c")d<. (29) 



Apparently, the only formal difference between 



(28) and (29 1 is the different measure. 



Due to the conformal equivalence, for a holo- 
morphic function / on the closed upper half plane 
Im z > 0, under the assumption that there is an 
a > for which |z"/(2)| is bounded, one has the 
Schwarz representation in the following form: 



/(.)^ir ^--^(-+^Q")d.. (30) 

^ ' TTl/ _ X-Z ^ ^ 



Note the formal equivalence to the spectral rep- 
resentation of a conventional Matsubara Green's 
function, Giz) = J dx ^ ^^^^^^±i^ . 

Similarly, the Poisson kernel for the closed up- 
per half plane is 



x'^ + 



(31) 



yielding the representation 



1 



u{x + iy)^- Pyix-t)f{t)dt, (32) 

with f e LP{R). 

The full Herglotz-Nevanlinna representation of 
arbitrary analytic functions with positive real part 
for the open upper half-plane read^^Hl 



/(^) =- 



z + i f dfi(x') 



Tri J (x' — i){x' — z) 



1 f d/i(.T') 



(33) 



TT 1 + X'- 



iaz + b. 



Here, ^ is given by the boundary-value distribu- 
tion of Re /: 

/i^Rebv/, (34) 
i.e. fj,{x) "=" Kef{x + iO+). The linear coefficient 

1 f d^{x') 



a = Rc/(i)-i| ^ 



„/2 



(35) 
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and the constant term b — Im/(i). 

For example, in the case of the function f{z) = 
i/vrz, /i is the Dirac measure fj,{x) = Rebv/ = 
6{x) and the coefficients a — b ^ 0. The case 
fj,{x) = —a-S{x) is not permitted by construction. 



At first glance, Eqs. (30 1 and (32 1 and the con- 



nection to the Herglotz-Nevanlinna representation 
( 33 1 seem to be rather straightforward applica- 
tions of the Cauchy integral formula. However, 
attempting the multidimensional generalization, 
we found that in our case, d > 1, the Cauchy- 
Bochner way of invoking JCc for a representation, 
equation (26), is not valid for the noninteracting 
Green's function GqIj^c: see Section 



HIE 4 



markably, as found in section |III E 5 



Re- 

taking as- 
sumption 2 from section |III D| in to account, we 



find that when transformed to the bicylinder, the 
Green's function Gj^c °iTic^ is closely related to a 
Caratheodory function, but the Cauchy-Bochner 
representation is invalid. 

These subtleties are apparently reflected by the 
central assumption (24 1 of Cauchy-Bochner rep- 
resentations. 



2. Functions with positive real or 



imaginary part in T 



Note that the representation (331 can be also 



used for bounded functions on the upper half 
plane. This can be seen by formally introduc- 
ing a shift in / which makes the real part of the 
function of consideration positive definite. 

While the signs of its real and imaginary parts 
will vary, the Green's function GLc is in fact a 
bounded function, assumption (12). This is why 



Vladimirov's integral representation for functions 
with positive imaginary part turned out to be 
applicable.^ Let us denote the set of holomorphic 
functions with positive imaginary parts on T^' by 
H^{T'~'). Due to the biholomorphic equivalence 
of the Green's function's wedges to the bicylin- 
der, one may think of H^{T'-^) as a generalization 
of the Garatheodory functions. Note that in the 
literature, sometimes functions with positive real 
and sometimes functions with positive imaginary 
parts are considered, resulting in marginal differ- 
ences in the equations. 



3. Vladimirov's Kernel Functions for T 

We will now study the generalization of 
the Herglotz-Nevanlinna representation to d- 
dimensional wedges ."^^ Vladimirov's approach 



may be found for positive real parts in Ref. [20] 
and for positive imaginary parts in Ref. 1131 

Because our original work^^I referred to Ref. \W\ 
we would like to switch to considering the class of 
functions with positive imaginary part, H^{T^), 
in the following. 

Let us first introduce Vladimirov's generaliza- 
tions of the Poisson and Schwarz kernels, using 
the Cauchy kernel ICc from Eq. ( 26 ) as a starting 
point. 

The (generalized) Poisson kernel for the wedge 
is defined bjll^l 



z = x + ly. 



(36) 



In case of the tubular octant, it is the product 
of usual Poisson kernels, however it is no longer 
simply proportional to the imaginary part of the 
Cauchy kernel. The (generalized) Schwarz kernel 



with respect to a point z*^"-* 
given by 



.(0) 



(0) g 



(0) . 2/Cc(z)/Cc(z(0)) 
For a measure /i(a;) we call 



(37) 



Pc[M{z) j M^)'Pc{z-x) (38) 
the Poisson integral with respect to /i. 

4- Vladimirov's Theorem 

We already pointed out that in the case d > 3 
two arbitrary different wedges, and , are 
usually not biholomorphically equivalent. In par- 
ticular, the wedge is not necessarily biholo- 
morphically equivalent to when d > 3. Hence, 
one may expect the structural similarity to the 
Caratheodory functions to break down more eas- 
ily in higher dimensions. This is the reason 
why Vladimirov's d-dimensional generalizatiorfi^ 
of Herglotz' theorem is stated in a comparably 
cryptic way which will simplify considerably in 
our case d = 2, for the reasons above: 

Theorem. (Vladimirov, 1978/79) The fol- 
lowing conditions for a function / e H+{T'-') 
are equivalent for a cone C C M*^ and ijl{x) := 
Im/(x + iOC): 

1. The Poisson integral Ppfd/i] is plurihar- 
monic in T*^; 
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2. the function Im/(z), z = x + iy £ , is 
represented by the Poisson formula 



lmf{z) = Pc[dfi]{z) + {a,y), 



(39) 



for some a £ C* , where C* is the dual cone 
of C; 

3. for all z^ € T*-^, under the assumption that 
C is regular, the Schwarz representation 



/(z)=i/ 5cU-t,z°-t)d^(i) 



(40) 



+ (a, 0) + 6 



holds, with b = 6(zO) = Re/(zO) - {a,x°)n 

Note that pluriharmonic functions are the natural 
multidimensional generalization of harmonic func- 
tions. A regular cone C in our context is a cone 
for which l/K-c is non-singular in T*-^. In the cases 
d= 1, 2, 3 all pointed cones are regulail^Sl. 

Using the equivalence of all T*-^ in the case 
d — 2, in Ref. [TU| we verified that the first state- 
ment of the theorem is true for G\rpc- This is so 
because it is known from the literature (see p. 134 
in Ref. 20) that the Poisson integral is plurihar- 
monic for any function H+IW^). Due to the bi- 
holomorphic equivalence of all T*-^ to in C^, 
the two integral representations provide exact re- 
lations for all holomorphic sheets of the interact- 
ing Green's function. A parametrization of the 
Green's function with respect to their Edge val- 
ues is gained by this. In our case, the validity 
of the representation has due to the biholomor- 
phic equivalences first been shown by Koranyi and 
Pukansky's work on the polycylinder.'^^ 

Note that the linear growth term a is zero for 
the Green's function, because it is bounded, as 



required by assumption ( 12 ) 



G. Application to the Green's function 

It turns out to be reasonable to specify a given 
cone domain arising from the branch cut structure 
by an angle 'd and an opening ratio r. See figure 
[TS] It is sufficient to consider the case •& = first, 
because relations for finite may be reconstructed 
from biholomorphic rotations, as explained in sec- 
tion UnEll 



1. Kernels for i? = 

For the case i? = 0, a computation of the ker- 
nels ICcr a and ^) has been provided in the 




opening 
ratio r 



FIG. 15: An arbitrary cone C with as vertex 
can be parametrized by an opening ratio r and an 
orientation angle z9. The rotation which maps 
Crfl to Cr^-d induces a biholomorphism between 
T'^'--° and T^'--'' . 



appendix of Ref. \W\ already, where the opening 
ratio r was named £ for technical reasons.'^ 
We used the definition 

Cr,a U {(a;i,a;2) e |x2 > OAxi = Axa} 

Xe(-r.r) 

(41) 

and computed the Cauchy and Poisson kernels via 



equations (26 1 and (36 1. The resulting Cauchy 
kernel is 

/Cc.,ofe) = -2r n ^ ' (42) 

and the resulting Poisson kernel is 



n 



ry2 - iJ-yi 



X±i (™2 - ^^XlY + (ry2 - m)^ ' 

(43) 

We have not used an explicit formula for the 



Schwarz kernel (37) yet, because the occurance of 
the reference point z*-*^^ appears to introduce addi- 
tional technical complications. The shape of equa- 



tion ( 43 1 is so simple because it can be computed 
from the tubular octant, whose Poisson kernel is 
the product of usual Poisson kernels. A simple 
real-valued 2x2 matrix acts as the biholomor- 
phism which converts the two representations. 



2. Operator Notation 

In order to put a stronger emphasize on the 
functional-analytic nature of the integral repre- 
sentations which interrelate edge and wedge val- 
ues of the Green's function, let us introduce an 
operator notation for the Poisson integral and also 
for the biholomorphic rotations. 
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Let us denote the set of all a-priori admit- 
ted Green's functions on T'^ by Gc- By "a- 
priori admitted" we mean those analytic functions 



for which axiom 2 from section 



HID 



holds (T*-^ has to comply with axiom 1). Fur- 
thermore, let us denote the corresponding space 
of boundary value distributions (edge functions) 
G\^c fe + iO^) hy£c- 

In order to focus on the Poisson kernel, we 
introduce the corresponding spaces of imaginary 

parts, Gc ^ and £q ^ ■ 

We denote by the operator Vc the linear map 



Vc : -> 



(44) 



ImG(x-t-iO'^) h-> dx^Vciz- x)lmG{x + iO'^'). 



Note that the Schwarz kernel (37) does not di 



rectly yield a comparable map from 8c to Gc: 
due to the occurrence of the z^^^ reference point. 

Furthermore, the rotation R^, which maps the 
cone Cr^o to the cone Cr,-d, induces a biholo- 
niorphic map : T'~''''-° — >■ T'-^^.^i (see section 



III E 2 ) . At this point we would again like to em- 



phasize that the biholomorphism does not con- 
nect the different branches of the Green's function 
on the wedges. It merely yields a counterpart of 
a given holomorphic branch on a biholomorphi- 
cally equivalent wedge which can be formally op- 
erated with. It is in that sense that is analogous 
to the concept of a conformal map. This biholo- 
morphism maps functions / e Gcr o to functions 
/ € GCr a • This can be similarly represented by 
the linear operator 



In practical computations, the function can be 
evaluated combining equation ( 43 ) and the rota- 
tion matrix 

„_i f cos-d — sinz?\ 

=(sm^ cos^ )• (49) 



4- Edge Properties of Go 

Since we essentially reduced the structure of the 
Green's function to the edge values of their holo- 
morphic branches, it seems worthwhile to investi- 
gate the edge structure of Go, and later also the 
perturbative structure of the theory in U, more 
carefully. See section [III I| for the deeper analysis. 

The edge limit of the bare Green's function (|8|, 
as a function of the edge orientation i9, is given by 



G. 



(edge) 




E 

a=±l 



G. 



(edge),, 



'i&^x), (50) 



where 



G, 



(edge) , a 




id;x) 



a{x^ - $)/2 - £d + irsgn^ 



and 



sgn^ := sgn ( cos ^ ^ sin d 



(51) 



(52) 



Apparently, the edge function only changes as a 
function of i? whenever cos iJit sin 19/2 crosses zero. 
This reflects the equivalence of all directions 



i^Co ^ Qc.,., (45) 
f{z) ^ f{R^\z)). (46) 

The operator TZ^ also naturally extends to a lin- 
ear map from o to a which we will denote 
by the same symbol TZ§. 



iO^ X 



sinz? 
cos-d 



0^ 



(53) 



when approaching the edge within a h olomor- 



phic branch T*^ , as discussed in section III E 2 



The edge function changes whenever •& crosses a 
branch cut, namely for the following singular ori- 
entations in the interval [0, 27r): 



3. Kernel functions at finite i9 

Consequently, for finite i9, the Poisson kernel 
operator of T^^-^) is 



(47) 



Equivalently, the Poisson kernel function of T'-^'-<^> 
is given by 



(48) 



(sing) 



^(sing) 
^(>5ing) 
^(sing) 



arctan2; 
TT — arctan2; 
TT + arctan2; 
2tt — arctan 2. 



(54) 



These are the angles corresponding to the four 
half-lines emerging from the origin in figure |12a[ 
The orientations -jj^'^'^s-' are also identical to the 
singular directions of assumption 3' in section 

inrDi 
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There is another subtle feature of the edge be- 
haviour of the bare Green's function. The real 
part 

ReG[,°''^^^(i9;x) = 

^{xu, - a{x^ - $)/2 - Ed) (55) 
at^i (a;c.-a(x^-$)/2-e,)2 + r2 

is completely i^- independent. As a consequence, 
for any branch of Go, the edge limit Re Go(a:+iO"-^) 
is identical. 

Another property is that, following the instruc- 
tions ( 10 1 to obtain the physical limit as far as the 



orientation of the limiting procedure is concerned, 
the function 



Ao(x) := --ImGo (x + i-0 
vr ^ 



(56) 



is positive definite: 



Mx) = E 



(57) 

In particular, the non-interacting spectral func- 
tion 



Aa{uj) = io(*,'^)- 



(58) 



adopted to the inference of A and therefore the 
spectral function. 

The procedure was found to work well in the 
equilibrium limit, $ = 0. However, entering the 
nonequilibrium regime, the ill-posedness of the 
inverse problem increased. Similar to the inter- 
twined geometric dependencies between the func- 
tion structures on edge and wedge coming to the 
surface in the appendix' uniqueness proof, a geo- 
metric dependency of the quality of Bayesian in- 
ference was found. 

Decreasing the parameter e provided a limit to 
a holomorphic function (leaving aside the vertex- 
correction branch cut) on the one hand, but on 
the other hand increased the ill-posedness of the 
inverse problem for a finite-$ spectral function. A 
discussion of how this is reflected by the structure 
of the Poisson kernel function may be found in 
RelM 

Apparently, the problem is very much related 
to restricting to the sheet T^^^" only taking 
G{iipm, it^n) data from the sheet into account, dis- 
carding the others. 

The only possible way to alleviate the increas- 
ing ill-posedness is to provide a link between the 
holomorphic branches of the Green's function, be- 
ing able to take into account data from not one 
but several wedges in order to perform the ana- 
lytic continuation (10). 



Again one can see in equation ( 57 1 that Aq does 



not decay to zero as a function of a; — >■ oo along the I. Bogolyubov's edge-of-the- wedge theorem 



singular directions (x^'^^ in axiom 3', section 



HIE 2). This is because the singular directions 



are an essential feature of the edges and lead to 
a discontinuity at oo when one compactifies the 



edge as shown in figure 14 



H. Bayesian Inference of Spectral Functions 

In Ref. [To] we used Vladimirov's integral rep- 
resentation in order to reconstruct a function A 



which was defined by equation ( 56 1 for the inter 



acting system. We chose a cone domain with ori- 
entation zero, T*-^^ ", and assumed the constrained 
Green's function G|^c^_o to be analytic for suf- 
ficiently small cone opening ratios e. This was 
justified, because the higher-order branch cuts of 
particle-hole character (see figure 12b I occur only 
in high-order terms in U. The vertex-correction 
type of branch cut pointed out in Ref. was ig- 
nored. 

Then the standard maximum entropy 
procedurd^ for inferring spectral functions 
from quantum Monte-Carlo data could be 



A candidate of such a link was provided by Bo- 
golyubov's famous edge-of-the-wedge theorem in 
the context of axiomatic quantum field theory. It 
considered the analytic continuation of Wightman 
functions^^ in order to establish certain dispersion 
relations. From a mathematical point of view, it 
also introduced a generalization of the very notion 
of analytic continuation.^^ 

There are several versions of the theorem. A 
simple version which captures the essential idea 
may be found in the book by Hormander on par- 
tial differential operators^^. It roughly considers 
two functions which are holomorphic on the 
tube cones r±'=', where G IS a convex open cone 
with vertex at zero. Consequently, the edges of 
T*-^ and T^*-^ are "infinitesimal neighbours". If 
the functions have the same boundary value dis- 
tributions, /+(x + iO<^) = /"(^-iO<^) =: /o, then 
/o is an analytic function, /o provides an analytic 
continuation of both, and /~ . 

In its more general formulations, the theorem 
actually demands the functions to be holo- 
morphic only locally at the edge and establishes 
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certain facts about the domain in which /o is an- 
alytic (global edge-of-the- wedge theorem). 

An extension to several cones whose edges meet 
in a single point is Martineau's theorem. As in 
Bogolyubov's theorem, locally, holomorphic func- 
tions may be found which constitute analytic con- 
tinuations of pairs of functions living on wedges. 
Again, the edge values of the considered set of 
functions have to be interrelated in a more or less 
direct way. 



IV. SYSTEMATIC EXTENSION OF THE 
CONTINUATION PROCEDURE BY USE 
OF EDGE RELATIONS 

In the previous section, we systematically an- 
alyzed the function-theoretical structure of the 
Green's function with regard to the two complex 
variables and z^p. The former comes along with 
the analytic continuation of the fermionic Mat- 
subara frequency ia;„ associated to the dynami- 
cal properties of the effective-equilibrium systems. 
The latter comes along with the analytic con- 
tinuation with respect to the Matsubara voltage. 
A fundamental property of the Green's function 
with regard to the two variables is the branch cut 

It separates the 



structure shown in figure 12b 



holomorphic sheets of the Green's function which 
live on wedges. Their edges meet in a branch 
point. For a holomorphic sheet, we were able to 
derive an integral representation of the Green's 
function with regard to two real variables, using 
kernel functions such the one in Eq. (43). By this. 



the functions values on the wedge are represented 
linearly by boundary values on its edge and vice 
versa. The physical limit ( 10 1 of the theory corre- 



sponds to approaching the branch point in figure 



12b along a certain direction. In order to use data 
from several wedges for physical results, it is thus 
necessary to find more or less explicit relations be- 
tween function values on edges of different wedges. 
The so-called edge-of-the-wedge theorem (section 
nil) provides some insight along this line. In or- 



der to construct an explicit functional-analytic ap- 
proach to the analytic continuation which would 
enable us to extend the numerical implementa- 
tion of the maximum entropy approach, it is only 
of indirect use, however. 

It is clear that any simple relation between 
edges of the different branches of the Green's func- 
tion provides a rather direct link between integral 
representations of the respective wedges. Based 
on a continuity approximation to function values 
at the branch point around which the edges are 
aligned, the present section derives the MaxEnt 



procedure which was used to infer the numerical 
results of section In] 



A. Continuous real part at branch point 

The relation is an exact identity of the bare 
Green's function. Namely ReGo(a; -f iO'') is iden- 
tical for any edge orientation d of the bare Green's 
function, see equation (55). It is easy to ver- 



ify that same is true for the second-order self- 



energy (151 and also for the functions which are 



parametrized by Han and Heary's original fitting 
approach in reference ilj 
We have 

ReGo(a; + iO'') = Y^cGo{x + iO^'), (59) 
ReS(2)(^ + iO'') ^ ReE(2)(^ + io^')^ (go) 

for all I?, t?' e [0, 27r). This structure is similar to 
the conventional Green's function causality rela- 
tion 

G{z*) = G{z)*. (61) 

There, we consequently have 

ReG'(a; + iO+) =RcG(cj-iO+). (62) 

However, in our case we only know for sure the 
symmetry 



G{z* , z* ) — G(z^p, Zi^y 



(63) 



From equation ( 63 ) only the edge relation 

ReG(x + iO''+'') = RcG(a: + iO''), d e [0,tt) (64) 

can be derived. I.e. conjug ate wedges , T-^ 
carry the same real parts of G(z) on their edges. 

B. Range of the continuity assumption 

We will now investigate to what extent the re- 
lations (59) and (60 1 also hold for higher-order 



contributions to the fully interacting Green's func- 
tion G(z), i.e. to what extent we a-priori expect 
the approximation 

ReG(^ + iO'') «ReG(x-HiO''=°) (65) 

to hold. It is insightful to study the algebraic 
properties of a conventional Green's function G{z) 
first. Subsequently, the two- variable function 
G{z) is discussed with respect to its formal struc- 
ture and regarding empirical findings from the 
continuous-time QMC simulation data. 



25 



a. Conventional Green's function. As a sim- 
ple examplejlet us consider the summation of the 
Dyson series^ 

G = Gn-GnJ:G (66) 
= Go — GqSGo + GqSGoSGo • • • . 

The entities G(z), Go{z), and satisfy the 



causahty relation (61). For the equation to hold, 



the product (and the sum) of two (61 l-satisfying 



quantities A{z), B{z) shall also satisfy (61). This 
is obviously the case, because 



Rc{AB){uj + iO+) = 



Re A{uj + iO+)Re B{uj + iO+) 
-lmA{uj + iO+)lmB{oj + iO+ 
RcA{uj - iO+)ReB(a; - iO+) 
-(-ImA(w-iO+)) 
(-ImB(a;-iO+)) 
Re {AB){uj ~iO+). 



Same can be shown for the imaginary part. It 
is crucial to recall that the mutual conjugation of 
imaginary parts of upper and lower functions has 
been used explicitly for closedness of (61) under 
multiplication. In other words: the set of func- 



tions with only the property ( 62 1 is not closed 
under multiplication. 

b. Two-variable Green's function. The 
closedness under multiplication is in general 
violated for functions with solely a continuous 
real part on the branch point. This is due to the 
fact that no statement about the imaginary part 
is made, and case of the causality relation (61), 



conjugateness of the imaginary parts is needed 
for closedness of the real part's continuity under 
multiplication. For instance, one can easily verify 
that Go{z) ■ Go{z) yields different real parts on 
the edges. Same can be shown for Go{z)~^. 



C. Structure of the residual term 

Nevertheless, as discussed in this paragraph, we 
are able to show that the continuity assumption 
( 65 ) is recovered for a certain energy range. Ad- 



ditionally, empirical findings for the structure of 
CT-QMC data, as discussed in appendix [B} par- 
tially support the assumption by observing con- 
tinuity relations between edge functions. Last 
but not least, the assumption is justified a pos- 
teriori for a rather large collection of wedges via 
the obtained numerical results (see section IV E 
for further discussion). It is found in section 
IV E that including upper (and lower) wedges, 
|w„| > \(pm/2\ (and |aj„| < — 1(/3,„/2|) into the con- 
sidered collection of wedges causes the numerical 



procedure to fail. Otherwise, it converges. This 
observation is coherent with a strong violation of 
the continuity assumption at the principal branch 
cuts ImZ(^ = ±lm Zip/2. 

It is insightful to study how the resulting dif- 
ference of two given edge functions is structured, 
namely to study the local residue 

r(^^^'){x) :=ReG(x + iO'')-ReG(a: + iO''') (67) 

for arbitrary values {), d' G [0, 27r) as a function of 

c. Angular structure. Obviously, due to 



Eq. (64) we have 



0, 



(68) 



for ah R'^'^^'^^ = if ?9 and d'{+n) belong to the 
same wedge. 

d. Structure due to continuity of imaginary- 
time data. Using the empirical fact that the 
continuous-time quantum Monte Carlo data 
I](i</j„i7 iw,i) are continuous as a function of ipm 
and iOn, we can derive certain continuity relations 
for R^^'^ \ They are provided in appendix [b} 

e. High-energy structure. Let us also con- 
sider the high-energy limit 



|a;|»max{r,|[/|,|$|,|erf|}. 



(69) 



For this, \G{x 
zero only if Xi^ 



section IIII E 2 



-I- iO'')| is significantly larger than 
w ±Xip/2, according to point 3', 
The and "— " cases imply a 
separation of energy scales. For both of these two 
energy scales, the closedness under multiplication 
is recovered. 

One can easily see the recovery of the multi- 
plicative structure in the high-energy limit by in- 
vestigating the bare Green's function (50). The 
and "— " energy scales are given by the a — — 
and a = -\- addends (51 1 in Eq. (50), respectively. 



Each of the addends satisfies the multiplication 
rule, because the absolute value of their imaginary 
part remains the same for all "i? at fixed x. Same is 
true for the sum of the corresponding J2i cti — 
addendends in the second-order self-energy (15). 
Therefore, we conclude that 



lim R^^'^"> = for all 



(70) 



in contrast to the limiting behaviour of G itself. 

/. Consequences for the analytic structure. 
Hence, although the assumption ( 65 ) is appar- 



ently only approximate, the error R^^ ^(a;) which 
is being done by assuming the relation is localized 
around in the z-space. It is remarkable that the 
intermediate-coupling numerical data presented 



in section |TT] appear to be rather precise in the 
low-energy region, although the violation terms 
are a-priori expected to be strong at low energies. 
From the a-priori perspective, the assumption 



( 65 1 gives a correct picture of how the wedges 
are related in the high-energy range. When as- 
suming the relation, additional low-energy degrees 
of freedom have to be introduced in order to re- 
obtain an exact continuation theory (cf. section 
[y] ) . Empirical data discussed in appendix [b] indi- 
cate that these degrees of freedom are comparably 
well-behaved. 



D. 



Functional-analytic consequences of the 
shared-real-part assumption 



We will see that the continuity assumption ( 65 1 



leads to a complete description of the entire func- 
tion G{z) on all wedges only as a function of the 
single edge ImG'(a; -I- 10''^°). This is extraordi- 
narily attractive from a numerical point of view, 
because by this, the number of degrees of freedom 
when doing the maximum entropy inference is 
not increased, but all imaginary-time theory data 
G{iipm,iu}n) may be taken into account, without 
any a-priori constraint. As in the single- wedge ap- 
proach of Ref. [ini the spectral function can still 
be directly extracted from the MaxEnt result. We 
will see that for functions which comply with the 
assumption, it in fact alleviates the ill-posedness 
of the inverse problem, as desired. 

g. Construction of the kernel. Starting from 
equation ( 65 ) we can derive a representation of the 
Green's function with respect to Im G(a; -I- 10''^°) 
in the following way. First, we introduce the 
Hilbert transform operator H as 



TtJ X2 - 



(71) 



Then we can use condition 3' (section HID I in or 



der to apply the Hilbert transform for computing 
real and imaginary part from each other on the 
boundary of certain H-isomorphic complex lines. 
One can show that 



ImGLc,^ = Q^.i, •ImG(x + iO''=°), (72) 

TT 

where we introduced the operator 

Qr,fl := Vr^-ffUiinn-^n. (73) 

Let us also introduce a symbol for the right part 
of the operator sequence. 



=) 



Im G{x + iQ°) 




FIG. 16: (color online) Action of the operator 

— Q^^'^^'^'^- It translates between functions liv- 
ing on edges of wedges with two different angular 
orientations. The orientation d is the one of the 
considered data wedge, and is the orientation of 
the physical limiting procedure ( 10 1. Initially act- 



ing on the physical edge function, the consecutive 
formal operations which comprise — 2^°'^^°^ either 
change the angular rotation or leave it invariant, 
as indicated by the respective arrows. 



The action of -2^°'^*^°^ on ImG(3: + 10''=") is 
depicted in figure [16] First, the Hilbert trans- 
form H with respect to the x^^ variable yields 
ReG(x 10^^="). Then, via equation ([65]) it is 
identified with ReG(a; -I- iO''). In order to obtain 
ImG(a; + iO'') one formally has to transform to 

in the do- 
The in- 



the biholomorphic equivalent of G|, 
main T'-^^-o via the operator TZ^^ = TZ- 
verse Hilbert transform —H yields the imaginary 
part of the edge value of the function T'^'--o C, 
z I-)- {TZ^^{G\j,c^_i))){z). Transforming the func- 
tion back to the edge of T*-^' using TZ^ yields the 
result ImG(x-hiO''). 

Using the Poisson kernel Vr,i), 



Eq. (48), the 



Green's function is obtained in the desired wedge 
T'-^^.^ . The entire procedure is contained in Qr,^- 

h. Feature of Q^J^'^^'^'^ : decoding branch cut ge- 
ometry from single edge function. The opera- 
tor 2^,°'^^°'' is a map of pure "edge" character. 
Therefore, it is worthwhile to study it separately. 
Q^cdgc) well-defined for any square-integrable 
function f{x), no matter which orientation -d is 
considered. 

Considering an edge function which is compat- 
ible wit h the Green's function properties 1,2,3' 



(section III D ) , Q^^'^' is not defined for the singu- 



lar orientations (54). For example, a straightfor- 



(74) ward calculation shows that applying Q^^'^^'^^ step 
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by step to Go{x+iO'^~°) yields exactly the formula 
(50), with the switching behaviour (52) whose 



value is undefined for the orientations (54). The 



missing square-integrability of the edge functions 
along these directions is the corresponding math- 
ematical reason. In particular, whenever Q^"'^^'''' 
crosses a singular orientation of lmGo{x + iO^^'^), 
it exactly generates the jump in Im Gq (a; + iO'' ) as 
a function of 



Consequently, assuming (65) is correct, both 



not only the holomorphic structure, but also the 
complete information about the entire branch 
cut structure, namely the exact geometry of the 
branch cuts, are encoded in the single edge 
function ImG(a; + iO°). Same is true for the 
{■d = 0)-edge limit of the second-order self-energy, 
Im S*^^^ (x + iO*'), due to equation (60). It is always 
square-integrable, except for the directions Xi^ = 



-3/2x 



X(^/2, x^ — Xip/'2, x^ — 3/2a^ 



namely for the geometry of the 2nd-order branch 
cuts. 

For practical computations, we find that an ex- 
ploitation of symmetries of 2^°*^®°' is mandatory. 
Those are translational invariance and scale in- 
variance, but no rotational invariance: For the 
translation operator 



{%'f)ix) -.^fix-x!), 
and for the homogenous scaling operator 
(Aa/)(x) := X^fiXx), A > 0, 

we have 

^ 0, 



^(cdgc) ^ 
^(cdgc) ^ 



= 0, 
^ 0. 



(75) 

(76) 

(77) 
(78) 
(79) 



The proof of these commutator relations is pro- 
vided in appendix [Cj Note that because direc- 
tional scaling implies a shear and therefore non- 
conserved angles in the shapes of ImG(a; + iO*'), 
it is no symmetry of the operator, in contrast to 
uniform scaling. 



E. Numerical implementation of Qr,^ 

The numerical implementation of the kernel 
Qr^^ is nontrivial. Assuming that G{x + iO°) is 
sufficiently smooth, we can represent it by super- 
imposing localized test functions which span the 
space of edge functions. 



i. Integral structure of the mapping. Qr,-d in- 
troduces a quadruple integral. The first two inte- 
grals are the two principal value integrals which 
come with the Hilbert transforms. The second 
ones are included by the Poisson kernel Vr,^- The 
integrations are formally very similar to a se- 
quence of convolutions A * {B * {C * e)), where 
e is an edge function. A crucial point is that due 
to the distributional nature of both, the principal 
values and the edge functions, the associativity 
rule cannot be expected to hold for these con- 
volutions (see section 4.2 in Ref. Principal 
value and our type of edge functions (functions 
with singular directions, Eqs. (54)) are no distri- 
butions with compact support. Therefore, it is 
impossible to simply contract some "inner inte- 
grals" within Qr,i? analytically in order to obtain 
a simple kernel function for Qr^§. The use of a set 
of test functions which spans the space of edge 
functions is mandatory. 

j. Construction of the test functions. The 
test functions would preferably be structured in 
a way which allows the quadruple integral in the 
operator Qr,i5 to be solved essentially analytically. 
Using the translation operators 7x , Eq. ( 75 ) , and 



scaling operators A^/^, Eq. (76), we define the 
functions 



with 



/(^) -in 



a— ±1 ^ -ri J 



1 



(80) 



(81) 



They turn out to be a good choice as test func- 
tions: First, we have the Dirac delta distribution 



8{x - X) 



(82) 



as a limit. Second, due to the symmetries (77) 



and ( 78 ) , the use of scaling and translation oper- 



ators yields - regarding the action of the integrals 
in Q^'^'^^"'' - the much more simple expression 



rather than (7'.,^(Qi"''^'VxAi/e/))(x) 
trix element of Q,r,§. Third, the simple pole 
structure of (IsTl) allows us to compute most of 



(83) 



as a ma- 



the integrals analytically. Note that the simpler 
looking symmetric Lorentzian function ^2 j^^i 
is, in fact, no good alternative to /, because the 
poles with respect to z^) or z^p contain square roots 
of z^ or z^j, respectively. Similar problems arise 
for localized Gaussians. 



28 



The directional arbitrariness x^^ ± aris- 



ing in Eq. (81) from choosing a product of 1- 



dimensional Lorentzians in (81) is still to be dis- 



cussed. For example, one could also have chosen it 
to be Xui^x^, adjusting the normalization factor 
from ^ to ^ in order to assert ( 82 1 . A concep- 
tional advantage of our choice of / is, however, 
that for any domain T'^'' '', for which G is holo- 
morphic, we have r < 2. Consequently, due to 



the pole structure of (81 ), / is holomorphic in the 
domain T'^^'° c T'^^-", whose edge is the starting 
point of the Qr^§ transform. Nevertheless, a cer- 
tain ambiguity remains which could be technically 
useful. 

k. Computation of the matrix elements. We 
found it feasible to calculate at least the first three 



integrals of the right-hand side of expression ( 83 ) 



analytically, using a computer algebra system. In 
order to compute the fourth integral, an adaptive 
numerical quadrature can be used. 

The result of the analytical integration of 
the first two integrals, namely (— is 
shown in figure [TT] for selected edge orientations. 
We find that {—Q^^'^^^^f){x) is a rational func- 
tion which changes continuously as a function of 
1?, in contrast to the transformation behaviour of 
ImGo(^ + iO°). Note that since Q^^'^^"^ is scale- 



invariant, Eq. (78), the transformation behaviour 
of the Dirac delta distribution is analogous to fig- 
ure [it] Consequently, the transformed delta dis- 
tribution on the edge M -I- iO'' is not a function 
but rather a distribution with a relatively com- 
plicated structure. Therefore, the limit 5 — in 



( 83 ) cannot be taken before the last two integrals 



from the Poisson kernel Vr,^ are computed. 
In the special case d = 7r/2, figure (17c), the 



asymptotic behaviour of the result decays oc j^y 
when a; — > cx), in contrast to the original test func- 
tion behaviour f{x) oc This is because the 
Hilbert transforms are taken with respect to mu- 
tually orthogonal directions in M? , here. 

The angles between and 7r/2 interpolate 
smoothly between the extremal cases of the un- 
perturbed well-localized /(x) at = and the 
long-range function at -d = tt/2. The solution 
at I? = TT is again strongly localized and equals 
~f{x). The behaviour in the interval [tt, 27r) is 
analogous due to symmetry reasons. 

/. Implementation. As mentioned above, the 



third integral of the operator sequence (83) can 



still be computed analytically. However, each in- 
tegration of the sequence adds additional poles to 
the resulting function, and more and more com- 
plex distinguishments have to be done in order to 
decide whether a pole is on the upper or lower 



half-plane and whether it contributes or not to 
certain residue sums. 

Because translational and scale invariance do 
not seem to be as useful concepts as applied to 
Vr^T}, not only the extra variable r appears in 
the computation of the remaining expressions, but 
also the shift X and the scale A of the test func- 
tion (80). For the third integral, one can still de- 



termine the poles and residues before doing the 
latter substitution with the computer algebra sys- 
tem, however. 

At present, very lengthy expressions result for 
the last integrand. As a consequence, the last in- 
tegral was evaluated numerically for each matrix 
element. The limit ^ — > can only be taken nu- 
merically, depending on the grid. An algebraic 
determination of the poles of this expression is 
cumbersome, because high-order polynomials ap- 
pear in the denominator of the resulting expres- 
sions. Nevertheless, numerical computations in- 
dicate that the limit (5 — yields well-defined 
functions after the fourth integration. Once an 
algebraic expression is found, the expression for 
the limit (5 — )■ would be more simple than the 
intermediate terms. As already stated in the be- 
ginning of this section, we compute the fourth 
integral with an adaptive numerical integration 
routine in practice. In the numerical MaxEnt im- 
plementation, one can adjust 5 as a function of x, 
denoted by 5x, depending on how well a specific 
region of the edge should be resolved. 

When defined according to the interacting 
branch cut geometry, Qr.i? will be simply called 
Q, in the following. Details on the numerical rep- 
resentation of Q are provided in appendix [Dj A 
detailed description of implementation and setup 
of the maximum entropy method using Q is given 
in appendix [E} 

m. Details of the MaxEnt procedure The 
quality of the results of the MaxEnt method us- 
ing the Q-mapping critically relies on a careful a- 
posteriori identification of the high-energy struc- 
ture of the physical edge function A(x^^ x^) along 
its singular directions. The information is in- 
corporated into the default model as described 
in appendix [Ej In brief, we first determine the 
most probable lateral width of the default model 
( E3 ) . In the same fashion, a second step optimizes 



the low-energy structure of the default model, by 
a posteriori determining the most probable low- 
energ y ba ndwidth, i.e. the quantity adef in equa- 
tion (E4) with the highest posterior probability. 



In equation (E4) we set the low-energy scale to 
i? sr. In tested examples, no strong depen- 
dence of the inferred results on i? was observed. 
In future applications however, in order to in- 
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(a) 1? = 



(b) i9 = 7r/4 



-Q 



{eda 



7W 



(c) 1? = 7r/2 



-Q' 



(d) = 37r/4 



-Q 



(edg 



7W 




(e) ^: 



FIG. 17: (color online) Transformation behaviour of the test function f{x) as a function of the edge-to- 
edge map —2^°'^^°'' for different values oft?. Function values are shown within the range [—5, 5] x [—5, 5]. 
Due to translational and scale invariance, it represents the edge-to-edge transformation behaviour of a 
Dirac delta function under the continuity assumption ( 65 ) . 



■ single analytic wedge (r=2, ^0) 
Q-approach 



FIG. 18: (color online) Comparison of as a 
function of the MaxEnt regularization parameter 
for single-wedge kernel 'Pr,{) and multi-wedge ker- 
nel Q at weak interaction U — 2T. (3 — SF^-'^, and 
e$ = F as a function of the regularization param- 
eter a. For the same input set, the single-wedge 
approach clearly fails to converge due to the pres- 
ence of higher-order branch cuts. 



crease accuracy, it may be advantegeous to also 
perform an optimization with respect to the pos- 
terior probability of R. 

Let us provide an example of how the method 
works for the weak-coupling case, i.e. for parame- 
ters U = 2T, 13 = 5F^i, e$ = F. The final resuhs 
for the spectral function were already presented in 
Fig. [3j It turns out that the applicability of the 
Q-approach is limited a posteriori by bad behav- 
ior of the inferred spectral functions to input data 
with w„ > |(/9„j/2| for n > 0. This condition cor- 
responds to not crossing the principal branch cut 
7 = ±1 in figure 12b when coming from the re- 
tarded Green's function (edge orientation = 0). 



Apart from this restriction, there appears to be 
no further problem with the approach. Conse- 
quently, at least for weak coupling, the central 
continuity assumption of the Q-approach is prac- 
tically solely violated with respect to the branch 
cuts 7 = ±1. The violation already occurs at 
very small values of the many-body interaction. 



U <2V. However it vanishes at [/ = 0, since the 
continuity assumption is exact for Gq. This ob- 
servation is compatible with the observed strong 
violation of the assumption within the Dyson se- 
ries which was reported in section |IVB| 

We therefore can use all Matsubara data of our 
Monte-Carlo simulation subject to ujn > \^Pm/2\. 
For the inverse temperature j3 ~ 5F^^, these 
extend from n ~ 1 to n = 8 for ujn with 
m = ±1,...,±5 for Lpm- As a first test, we 
show in Fig. [18] the performance of the Max- 
Ent method for both, the single-wedg^^^l and 
the multiple-wedge approach for the given data 
set. Because data from the comparably widely 
opened wedge a;„ > |(/9,„/2| are used in the single- 
wedge approach, it implicitly assumes the inter- 
acting Green's function to be analytic for Im > 
|Im2<p/2|. Apparently, this wrong assumption 
makes it impossible to obtain a reasonable fit 
with a positive definite A{x). Consequently, the 
value of the procedure does not drop below 
10^ • A^data, and the MaxEnt fails to converge. In 
sharp contrast, values of x^/^data ~ 1 may be 
reached with the MaxEnt with respect to Q. Also, 
controls such as the MaxEnt error rescaling merit 
do not indicate the presence of any abnormalities. 

Thus, for the Q-mapping, a well-behaved Max- 
Ent solution is obtained. As further discussed 
in appendix [Ej the quality of the solution very 
much relies on appropriately including the prior 
knowledge on singular directions of the Green's 
functions in into the default model of the edge 
function. For this purpose, within a set of smooth 
default models with the correct singular direc- 
tions as a; —> CXI, a most probable one is identi- 
fied within the Bayesian framework of the MaxEnt 
method. The thus identified default model for the 
edge function is displayed in figure 19a Using 



this default model, the well-behaved edge func- 
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FIG. 19: (color online) Application of the MaxEnt procedure for the Q-mapping to the nonequilibrium 
weak-coupling case U — 2r, /3 — 5r"^,e<i> — F, with CT-QMC data as input. The default model has 

been identified via its maximal posterior probability. 



tion A{xip,x^) shown in panel 19b 



An overall moderate sharpening of t 



is obtained, 
le edge func- 



tion along the cross-like structure is observed as a 
result of this final step of the Bayesian inference 
procedure. 

With such appropriately optimized default 
models, the results presented in section [ll] were 
obtained from weak to intermediate coupling 
strengths. Throughout, the only major data 
range constraint, w„ > |iy5m/2| was found, which 
prohibits crossing the principal branch cut due 
to violations of the continuity assumption. On 
occasion, for stronger correlation strengths, val- 
ues with small a;„ had to be discarded in order 
to obtain a converging MaxEnt solution, i.e. a 
solution which meets the continuity assumption 
constraints. At the comparably small inverse 
temperature /? — 5F^^ used, calculations require 
only moderate computer resources, mainly due to 
the comparably small QMC data space of approx- 
imately 50 imaginary-time-theory data points. In 
general, the amount of data will grow quadrati- 
cally as a function of inverse temperature, due to 
the simultaneous presence of Matsubara voltage 
and Matsubara frequency. Additionally, at low 
temperatures, sharp features in the spectral func- 
tion and hence the edge function A will have to be 
resolved, requiring an enhanced grid refinement. 
Altogether, matrix sizes in the MaxEnt will in- 
crease substantially when the temperature is de- 
creased. In particular, the computational effort 
for the generation of an appropriate kernel ma- 
trix (cf. appendix [d| grows dramatically and the 
memory consumption of the MaxEnt itself poses 
a limitation at lower temperatures at the present 
stage of code development. Additionally, it is 
well-known that the resolution of low-temperature 
features with the MaxEnt method requires a care- 



ful Bayesian analysis based on higher-temperature 
data, i.e. an "annealing procedure" , involving a 
sequence of QMC plus MaxEnt runs for a reason- 
ably fine temperature gricP^^^H^ 



V. PERSPECTIVE: UNBIASED 
Q- APPROACH 

From a mathematical point of view, the under- 
lying continuity assumption of the Q-approach is 
only approximate, because in higher orders of per- 
turbation theory, terms which do not character- 
ize the full collection of wedges, but rather just 
isolated wedges or subcoUections of wedges, are 
generated. These terms are manifested in discon- 
tinuities of the real part of the Green's function 
at the branch point Imz — 0. 

In order to extend the Q-approach to the full 
nonequilibrium Kondo regime U > 27rF, e$ 
Tk, ~ Tk, one has to take these contribu- 
tions into account. This requires the considera- 
tion of the full analytic structure of the theory, 
i.e. the full set of edge functions. 

As a consequence, extra terms have to be added 
to the representation of G(z^, z^^) within the Max- 
Ent procedure. Obvious candidates for such de- 
grees of freedom are the residual imaginary parts 
of edge functions 

R,,{x) ■.= lmG{x + iO-^-)-TT-{Q'^y'A){x), (84) 

for the edge of the n-th wedge with orientation 
i?„. Because the Q-mapping is exact at high en- 
ergies \\x\\, the terms Rn{x) are essentially local- 
ized within a finite radius around 0. This range 
is expected to be of the order of magnitude of the 
energy scales F, U, e^, and e$. 
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Regarding the inverse problem, for data in the 
n-th wedge, one has the exact representation 



Im G{i(pjn , iw„) = (Vr,, ,il„ Rn){i(Pm , it^n) 

+ vr- (QA)(i(p„,ia;„), 



(85) 



where r„ is the opening ratio of the respective 
data wedg^^. The MaxEnt procedure must de- 
termine Rn and A simuhaneously. Practically, 
the terms 7?„ would act as "valves" for the con- 
ceptual imperfection of the Q-mapping within the 
Bayesian information flow. 

It is an interesting question if the formally in- 
finitely many two-dimensional variable vectors in 
practice lead to a dramatic increase in the fit space 
or not. Due to locality of the terms Rn{x), the ef- 
fort is probably less than for the A function which 
itself encodes many aspects of the analytic struc- 
ture. Furthermore, the rather large Poisson kernel 
matrix elements at low energies will possibly lead 
to a comparably good MaxEnt performance in the 
determination of Rnix), as long the opening ratio 
?9„ of the n-th wedge is comparably large. 

Because the functions i?„(a;) cannot be ex- 
pected to be positive, it is necessary to introduce a 
shift to a positive function, such as for the spectral 
functions of the static observables in paper 1. The 
terms i?n(x) are presumably most dominant for 
wedges next to the noninteracting Green's func- 
tion's branch cuts. A very careful Bayesian anal- 
ysis, including an appropriate set of choosable de- 
fault models constructed from a-priori informa- 
tion, is probably required for a successful applica- 
tion of the exact approach (85). It is also possi- 



ble that the perturbative structure of the theory 
reorganizes terms i?„ in subcoUections of wedges 
which result in a more moderate MaxEnt problem 



than equation (85). In particular, the branch cut 



at Im z„ = probably leads to a nonzero limit 



lim„ 



, R„ , where the limit n 



oo shall con- 



sider a sequence of wedges with '&n^oo = or 



VI. SUMMARY 

We systematically studied the mathematical 
structure of the dot-level Green's function and 
the Bayesian inference of non-equilibrium spectral 
functions and transport properties from effective- 
equilibrium quantum Monte-Carlo data within 
the Matsubara-voltagc theory. Furthermore, a 
continuity assumption on the analytic structure 
was introduced which strongly improved the nu- 
merics of the MaxEnt approach of an earlier 
publication.'^ 



Formal parts of the paper introduced the es- 
sential concepts of the function theory of sev- 
eral complex variables and connected to the re- 
spective mathematical literature. Using insights 
from perturbation theory, the Green's function 
was characterized axiomatically with regard to 
its function-theoretical structure. As the funda- 
mental domains of holomorphy, so-called wedges 
(tubular cones) emerged. The Green's function 
is composed of sheets which are holomorphic on 
the wedges enclosed by branch cuts. Within each 
wedge, the Matsubara data of the Green's func- 
tion uniquely map to a real-time limit on the so- 
called edge of the wedge. For this purpose, an 
explicit integral representation was constructed. 
However, depending on the considered wedge, the 
edge structure does not necessarily have a direct 
physical interpretation. 

In an earlier publication,!^^ we had been un- 
able to compute reliable non- equilibrium spectral 
functions from a MaxEnt procedure based on inte- 
gral representations within wedges, due to rather 
strong assumptions on the analytic structure and 
weak assumptions on the physical structure. The 
assumptions had limited us to a single wedge with 
a rather small opening ratio, on which the Green's 
function is not strictly analytic but which directly 
includes the physical limit procedure on its edge. 
While as compared to the present work the nu- 
merical effort of the MaxEnt procedure was rather 
low, due to the wedge structure and simple ker- 
nel structure, we had not been able to consider 
most available quantum Monte-Carlo data within 
the MaxEnt procedure. The hereby implied loss 
of information from available simulation data had 
not affected the equilibrium spectra but the non- 
equilibrium spectra, due to the kernel structure. 

In order to overcome these previous limitations, 
we introduced a continuity assumption to the real- 
time structure of the Green's function, i.e. its 
structure at the branch point around which the 
edges of the wedges associated to the branches of 
the Green's function are aligned. The assumption 
includes structures generated by the earlier fit ap- 
proach introduced in Ref. [U which was motivated 
by perturbation theory. Mathematically, the as- 
sumption lead to a uniform description of data 
from all branches of the Green's function and gave 
rise to a linear operator Q which, while hard to 
implement, enhanced the MaxEnt procedure to a 
larger set of quantum Monte-Carlo data. 

We found that the continuity assumption ap- 
pears to be valid for a very broad range of data, 
up to intermediate coupling strengths, eventually 
yielding reasonable non-equilibrium MaxEnt re- 
sults for spectral function and transport proper- 



32 



ties, which are dramatically improved as com- 
pared to the results of Ref. [THl However, as 
the nonequilibrium Kondo regime is approached, 
we expect the continuity assumption to break 
down eventually. For this parameter regime, the 
method could be extended along the line discussed 
in section IVl 
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Appendix A: Uniqueness of the Analytic 
Continuation of Dynamical Quantities 

In the following, we would like to show that 
the continuation of Matsubara data G{i(pm,i^n) 
to the multisheeted holomorphic function G{z) is 
unique, i.e. we will prove assumption 3 in section 
HID relation (13). We will derive the unique- 
ness using the axiomatic statements 1, 2, and 3' 



of section |IIID| Since the proof will involve some 
elementary geometry, it will be accompanied by 
several sketches. 

We may focus our attention to a single wedge 
T*^ which is defined by subsequent branch cuts 
from Eqs. (111. The data G'(i(/3m, icj„) which 



are located in the wedge are our starting point, 
{ipjn,0Jn)'^ G C. Without loss of generality we 
can assume that we have entire lines of data, 
G{iipjnT^'-^i), '^i G I^, because arbitrary continu- 
ous imaginary uji may be computed by Fourier 
transform in the (^mth effective equilibrium the- 
ory, having again ((/3m,a;/)^ G C. Let us denote 
the effective equilibrium data range by 

Eq {i((p™, W7)^|to e Z, e M} n . (Al) 

These lines of known data of the unknown func- 
tion G{z) in the wedge T*^ are depicted in Fig. 20 
They constitute one-dimensional lines in the four- 
dimensional wedge T*^ for which the function 
G\rpc shall be reconstructed. 



FIG. 20: (color online) The wedge to be con- 
sidered. The dash-dotted lines denote the data 
yielded by imaginary-time theory. 



We will, step by step, prove the uniqueness 
of the continuation of the imaginary-time data 
by applying biholomorphic maps and the identity 
theorem of complex analysis. The central idea will 
be to extend larger and larger subsets for which a 
unique continuation is obtained. 



1. Reconstruction of edge values using 
complex lines which are isomorphic to EI 

Due to assumptions 1 and 2, we found that the 
Green's function G\rpc may be reconstructed from 
their edge values, using Eqs. ( 39 ) and ( 40 ) . There- 



fore, it suffices to show that we can reconstruct all 
edge values of the function G from the data Gj^^. 

We will first show that one may reconstruct a 
certain set of single lines through zero on the edge. 
Each of these lines is defined by an angle -i^. All 
function values on the line may be reconstructed 
if the angle d is contained by the cone G. 

The proof of the latter statement is the follow- 
ing. Consider a single line in the cone C, given by 
the angle d. A biholomorphic rotation in the sense 
of section [TlI E 2| can then be applied in such a way 
that the line is horizontal and may after complex- 
ification be interpreted as the upper half plane EI 



of C, see Fig. 21 The real line is then associated 
to a horizontal line on the edge of R-^T^ . 

The biholomorphic equivalent to the yet un- 
known function is now G{z) = G{RZI -z). Note 
that the line (0,iA)-'" in R^^T^ , A > 0, contains 
infinitely many known values of G{z). These are 
denoted by the crosses in Fig. [21] Extending to 
the upper half plane (0,IHI), one may apply the 
identity theorem of complex analysis for recon- 
structing G(z) on the whole plane (0,11). In par- 
ticular, the boundary values G{0, M -t- iO+) are re- 
covered. This proves the statement of this sub- 
section. 
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FIG. 21: (color online) A single line contained by 
the cone C may after biholomorphic rotation 
and subsequent complexification be interpreted as 
the upper half plane H of C. 



reconstruc- 
ted area_ 



reconstructed , 
area 




FIG. 22: 



Edge^^r"^"^ 



Al 



Uniquely reconstructed range 

of G{z)\j,c on the edge of 
following from the partial argument of subsection 
The wiggly lines in the boundary mean that 
the area extends to infinity. 



Matsubara Green's functions. 

One may also think of G'(C2:^°-'), C G lEI as 
a meromorphic function of C G C, because it 
may due to boundedness and 1/C asymptotics be 
approximated arbitrarily well by a meromorphic 
function, such as in an infinite Fade expansion. 
Since meromorphic functions on C are holomor- 
phic on the Riemannian sphere, the identity the- 
orem holds at the accumulation point oo. 



a. Identity Theorem 

Let us comment on the satisfaction of the as- 
sumptions of the identity theorem. Since (0, oo) is 
the accumulation point of the known data points 
on (0, H), i.e. of the "series of crosses in Fig. [21]', 
given by (0,]HI) n R-^Eq, we have to show that 
G(0, 1/z) may be extended to an analytic func- 
tion at z = 0. Once this is possible, the function 
is uniquely determined by the set of function val- 
ues. 

Combining assumptions 1 and 3', we know that 
G{Cs±°^) (in the sense of assumption 3') behaves 
like a conventional Green's function, because the 
singular case coincides with a branch cut which is 
by construction not contained by the wedge. Due 
to this rapid decay one may extend G to the lower 
half plane such that G{0,z*) = G{0,z)* and is 
holomorphic at z = oo. This can be done explic- 
itly using a spectral representation with respect 
to the boundary values of ImG(0,z) on the real 
axis. The spectral representation exists due to the 
1/z asymptotics which lets the line integral con- 
tribution vanish on the infinitely large semicircle 
attached to the real axis. Note that this construc- 
tion is also compatible with the symmetry relation 
G(-i(/3„j,-iw„) = G(i(/3,„,ii:j„)*. 

As a consequence, the identity theorem is appli- 
cable for G(0, z) at z = oo such as it is for regular 



6. Resulting reconstruction of edge values 

Sweeping through all possible angles which 
are contained by the cone G, the uniquely recon- 
structed edge behaviour of G\rpc is given by the 



area depicted in Fig. 22 



It is obviously given by 



Edge^™°) :=GU(-G). 



(A2) 



Extending the unique range to the entire 
edge 



In order to show that the function values of 
G\rpc are also uniquely defined by G\^^^ for the 

complement of Edge^c'^""'"'', the argument has to 
be extended in a similar way. The trick is to con- 
sider yet another set of H-isomorphic subspaces 
and then apply the argument of the last section 
to a larger set of data. 



a. Extending the known data range within the wedge 



The first step is depicted in Fig. 23 In contrast 



to before, we consider a constant angle iJo and 
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FIG. 24: (color online) Reconstructing the 
Green's function on the complete edge. The re- 

(see Fig. 



constructed area Edge^c'^°'^'°^ 



22 



and 

Eq. ( A2 ) ) may be extended by afhne translorma- 
tions along the i?^(,(0, 1)-^ direction using the in- 
formation from the set Ei. 



FIG. 23: (color online) Enhancing the formal 
holomorphic reconstruction within the wedge. 



various lines which start at different points on the 
boundary of the cone with the orientation i?o. 

After biholomorphic rotation to the wedge 
R-^gT'~^ we can again complexify the lines 

f.(o, := {iyW +i(0,A)^,A > 0}; y^'^ G dC (A3) 



to 



L^(o) := iy(") + 



(A4) 



The isomorphy of Ly(o) to EI and assumption 3' 
again enable us to apply the identity theorem to 
the crossed data in Fig. [23j namely to the infinite 
sequence (G o i?^)|(^_^g^)^£_^^^ . 

By this, the transformed Green's function G = 
G o R^g is reconstructed for all points of the set 



D :— Ly(o) — {R-^gT 



(A5) 



For simplicity, we may now just look at a subset 
of _D, namely 



El :— R-^ Eq 



(A6) 



It enables us to see that the values of G are now 
known on the set 



El :— i?o + Rda 



(A7) 



h. Full Reconstruction of the Edge 

With the information from Ei one may reinter- 
pret the procedure associated with Fig. [21] and de- 
scribed in section lA II The dashed lines of known 
data now contain an additional real dimension 
along the direction i?^(,(0, 1)-^. 

We can use each point A • i?ij(,(0, 1)-^ (A G M) of 
this new degree of freedom as an offset of the lines 



used in section (All and reapply the entire pro- 
cedure. Using the resulting afRne subspaces, the 
Green's function may be reconstructed on further 
regions of the edge which are afRne to the one in 



Fig. 22 namely 



Edge^^r"'^^ :=Edge(;;r"'"^ 



A-i?^o(0,l)^ 



(A8) 



Applying the argument to all A € M 
reconstructs the entire edge and hence 
the entire Green's function G|yc, because 

ASK Edge^c = Edge^^c . 

This "affine procedure" is sketched in Fig. [24] 



Appendix B: Empirical properties of the 
residual term 



The empirical observation that the quantum 
Monte Carlo data are continuous as a function 
of iuJn, as one crosses higher-order branch cuts, 
yields the following structure. Let be the ori- 
entation of the corresponding branch cut. Then, 
due to the observed continuity, we have 

G'(ipsin(i9o - (5),i/ocos(i?o - 1^)) = 
= G(i/9sin(i?o + (5), iyocos(i9o -I- (5)), 
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for any p > 0. 

Using the identity theorem along the directions 
i?o i 6, one find that the relation 

^(i>o-Mo+5)(^gin(^Q)^^cos(i9o)) = (B2) 

holds for any p > 0. Consequently, the continuity 
in Matsubara space induces a continuity relation 
in the edge space. Similarly, one could construct 
relations for the derivatives of G which help con- 
strain 



lated by the operator 73s:, where we set X = 
{X, Y)'^ . We also write x :— {x, y)'^ . 

As a first step, we apply the Hilbert transform. 



{HTxA){x) 



TTj y -y 

1 r_, A{x - X, y) 
tt/ {y-Y)-y' 



(CI) 



Appendix C: Commutator relations of Q^^"^^"^^ 



We will derive the commutator relations ( 77 ) 



(78), and (791 



Using the short-hand notation c cosi^, s := 
sin??, we then have 



Translational invariance 



(edge) 



Let us consider the action of Q}^ 

on a function A{x, y) which is trans- 



1^,^ 7x/ X 1^,- Alcx + sy — X,y) 
{n^^nTxA){x) = -Tfdy , ^ ^ ^ 

— n j {-sx + cy -Y) - y 

(C2) 

In the next two steps one obtains 




~ Ajjcx - sy)c + sy- X, y) 
-{cx — sy)s + cy — Y —y 



(C3) 



.{eg 



By substituting A 



and then applying Tx from the left, one finds 
I J 1 



dy 



A{{c{x - A) - s{y - Y))c + sy, y) 
s{x - A) + c{y -Y)-yJ -{c{x - A) - s{y - Y))s + cy - y 



(C4) 



In order to verify that the expressions ( |C3| and ([C4| are, in fact, equal, we consider the following 
system of linear equations: 



sx + cy — y ~ 
{cx — sy)c + sy ~ X — 
-{cx — sy)s + cy — Y = 



sxq + cyQ-y 
{cxq - syo)c + sy* 
-{cxo - syo)s + y*c, 



in matrix form: 




-1 s 






(y* 


s c^ 


— SC 




Xq 


c —sc 








= :M 







(C5) 
(C6) 
(C7) 



(C8) 



The equations correspond to the idea of substituting the terms in equation ( C3 ) as denoted there in 
such a way that y* is the new integration variable, J dy ^ k J dy* , where k is some regular prefactor 
from the integral transformation, and such that the external variables x and y are replaced by xq and 
yo- We will see that the resulting form will exactly be (|C4[): 
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The system (C5|-(C7) can be solved by inversion of the matrix M: One has 



M- 




(C9) 



detM = — 1, and consequently the well-defined solution 

— X — X, 
2/0 = y-Y- 



(CIO) 
(Cll) 
(C12) 



The integral transformation constant k—\, and performing the substitutions (C5|-(C7| in (C3| yields 



(C4), when y* is again renamed y. 



Therefore. 



0. 



2. Scale invariance 



Similarly, we show the scale invariance ( 78 ) . We have 



27'dy ~ 

TT ./ sx + cy — y 



dy 



. A{\[{cx - sy)c + sy], Ay) 



— {cx — sy)s + cy — y 
A{{cXx ~ sXy)c + sy,y) 
sx + cy — y/X j ~{cx ~ sy)s + cy/X — y/X 

d^ ^ 4dy ^^^^^"^ ~ ■'^^^'^ ^ 



sXx + cXy — yj —{cXx — sXy)s + cy — y 



(A.Q^f^'=)i)(x). 



(C13) 

(C14) 

(C15) 
(C16) 



3. Absence of rotational invariance 



We provide a simple example for which [Q'"^'^^°\tZ§'] 7^ 0. 



We consider the bare Aq, equation (57). Setting § = -d' — ti /2, we find 



and 



(7^V2Q^/T^)fe) = E 



^f^^ (x^ - a{-x^ - ^)/2 - saY + F^ ' 



(C17) 
(C18) 



On the other hand, 



{n^/2AQ){x) - E 



Fq/tt 



^ (x^ - - $)/2 - e<i)2 + F2 



and 



(Q(<^;f)7^./,^)(3:) - E 



-aFa/TT 



^ (x^ - a(-x^ - $)/2 - e<i)2 + F2 ■ 



(C19) 
(C20) 



Hence, Qt'/t°^7^,/2^ = -n^/iOi^^ff Ao, and therefore [Q^T^f \ 7^,/2] ^ 0. 
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Appendix D: Numerical representation of the multi-wedge map Q 

In this appendix, a recipe for the numerical computation of the quadruple integral Q is given. The 



application to the test function (80), fx.s^ for a function value at (i(y9m,i(u;„) is to be computed. The 
first three integrals can be computed analytically by use of a computer algebra system. A numerical 
quadrature method can be used for the approximation of the remaining integral. 

1. Analytic computation of first three integrals 

Using the translational and scale invariance of the edge-to-edge contribution Q,Jj''^^°'' , only the action 
of Q':^'^^'^^ on our test function 

/(^' ^ ^ _ ^/2)2 + l' {y + .t/2)2 + 1 



has to be computed, yielding the results shown in figure 17 For brevity we set y = and x ~ 



The first principal integral can be eliminated by straightforward application of the residue theorem: 



mf){x,y) ^Tf dy-^ f{x,y) 

J y~y , . 

{x'~12~Af)y ^ ' 

^ (16 y4 + 32 y2 _ 8 y2^2 + 16 + 8 a:2 + a;^) {x^ + 4) ' 

As a next step, introducing the short-hand notation s — sini^ and c — cos-iJ, imposing the rotation 
operator 72.^^, one obtains 

{n:g^nf){k, l) = subs {x^kc + is, y^lc- ks: {Hf){x, y)). (D3) 

Here, "subs" denotes the operation of a variable substitution. In order to apply the second Hilbert 
transform, it is necessary to determine the poles of the corresponding integrand 



(D4) 



g(k,y;l) :^ .{n^'nf){k,l) 

{k^c^ + 10 kds + Ps^-l2-4: fc^ - 4 k\s^) (Ic - ks) 
~ 7r3 {k'^c^ -I- 2 kds + Ps"^ + 4) 

• (4 Pc^ - 6 kds + 4 k^s^ +4 + 4 kc^l- 

-Ak^cs + 4fsc ~ Us^k + k^c^ + fs^y^- 

with respect to the integration variable I. One finds that the function has the following seven poles in 
the complex plane: 

h ^, (D5) 

-fcc + 2fcs±2i , , 

h.z ■■= , P6 

s 

-fcc + 2fcs±2i , , 

>?B±^, ,08) 

—s + 2c 

The Hilbert transform can now be evaluated through the residue sum 

7 

{WZ^^Hf){k,§) = [2^i-Res,=z„g(fc,^;/)] •0(ImZ„). (D9) 
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The Heaviside function d{x) ensures that only poles from the upper half plane are taken into account 
for the evaluation of the contour integral which corresponds to the principal value integral. An explicit 



evaluation of the residue sum ( D9 ) can be accomplished with a computer algebra system. A dramatic 
increase in complexity is coming along with the constraints Im^„ > which depend on the wedge 
orientation angle In fact, six separate cases emerge as a function of 'd. They can be parametrized by 
three overlapping different cases, namely (A) terms which are proportional to sgn(s), (B) proportional 
to sgn(2c + s), and (C) terms which are proportional to sgn(2c— s), discriminating between the different 
signs of Im^„. The term 

(7^^H7^^l•H/)(a;, y) = subs {k^xc^ ys, y-^yc^ xs; {nn:g^nf)ik, §)) (DIO) 

is then best reorganized into rational functions as coefficients of the sign functions, 

(Q^;''°V)(a:,y) = {n^mZ7/nf){x,y) =A{x,y) ■ sgn(s) + 

+ B{x,y)-sgn{2c + s)+ (Dll) 
+ C{x, y) ■ sgn(2c — s). 

This explicit split is necessary for the study of the interplay of the rational functions A{x, y), B(x, y), 
and C{x,y) in a computer algebra system. 

a. Rational coefficients of the transformed edge test function 
The rational functions 



A{x,y) = 


-^cnum('^5 y) 
-^dcnom ('^; 


(D12) 


B{x,y) = 


-^cnuin('^i y) 
-^denom('^i y) 


(D13) 


C{x,y) = 


^enuin (•^7 y) 
C'dcnom(-^; y) 


(D14) 



have the following polynomials as enumerators and denominators: 

/ N nf f 2 , 2 41 \ 4 13 /12 2 12 28 2\ -1 

^e„um(a;,y) =-8 I l-l + x^ + l/4y^-—xyjc*^— ( ^3 + ^ + 13 ~ 13 ^ ) '^'^ 

+ ^1 + y xy- l/2y2^ + l/16s + -Axy-Uy'^) c (D15) 
-l/l6y{-Ay + x)^, 

Adcnom(a;, y) ^TT^ {x^ +A-Axy + Ay^) ■ K{x, y), (D16) 



Bcnumix, y) = - 8 f (^-l + x^ + 1/4 + 11 xy^ c'^+ 



IS f 12 2 12 28 ,\ . / 13 ,„ 2A 2 

+ ^ . (- ^ + + ^ - ^ j + (l - y - 1/2 y^) (D17) 

- l/16s (-4 + a;2 +4xy- 12 j/^) c + l/16y (4y + x) ^ , 

-^dcnom (*^; y ) — 

(x^ + 4xy + 4y^ +A)tt^ ■ K{x,y), (D18) 
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Ce„um(x, 2/) = - 3 (-1/3 _ 3 + x^) yxc" - s ((-3 x2 + 4) „ 5 ^ 4 + x"") 
+ 3 (-10/3 + - 2/3 1/^) j/xc^ 
+ ((-3a;2 +4) ~ + 4) sc + yx {y'^ + l) , 



4^ ^3 



and 



C'donom(a:, y) = {x^ + 4) TT^ • K{x, y) . 
For the denominators, we introduced the shared polynomial 

K{x,y) = (2/"* + (-3-6a;2) 2/^ + 3 a;^ - 4 + x'') c"* 

- 4 (3/2 + x^ ~ y^) xsyc^ + (-2 ?/'' + (2 + 6 x^) t/^ + 4 + x^) 

- 4 (1/2 + 2/2) XS2/C + 2/ + y''- 

I 



(D19) 



(D20) 



(D21) 



b. Composition of the rational coefficients 



Contraction with the Poisson kernel 



In contrast to the actual edge functions 
Q(cdgc)j^ the coefficients A{x,y), B{x,y), and 
C{x,y) are comparably ill-behaved. For example, 
one obtains 



As a next step, one of the integrals introduced 
by Vladimirov's Poisson kernel (47) will be eval- 
uated analytically. For this, the pole structure 
of both, the Poisson kernel (48), and the edge- 



B{x,y)\^^o 



x + y 



TT^ (x2 + 4) (x2 + Axy + Ay^ + A) 

(D22) 

Apparently, the function diverges for x — > 0. The 
functions A{x,y) and C{x,y) are similarly struc- 
tured. However, the actual edge functions, as dis- 
played in figure [TTI have no such singularities, but 
do rather represent smooth deformations of the 
test function when the wedge orientation angle ■)? 
is tuned. As a consequence, the real singularities 
are cancelled as the rational functions are added 
up in equation (Dll). For a consistent further 



evaluation of the action of Q on the test function 
it is thus necessary to study each of the full com- 
binations (Dll) separately. There are six possible 



combinations, namely the sectors 



(D23) 



< "d < arctan2, 

< I? < vr- arctan2,(D24) 
(D25) 



(a) 

(6) arctan2 _ _ 

(c) vr — arctan 2 < < vr 

(d) TT < I? < vr + arctan 2, (D26) 

(e) TT -I- arctan 2 < i9 < 2t: - arctan 2(D27) 
(/) 2tt - arctan 2 < i9 < 27r. ''r>'^«^ 



(D28) 



For example, the expression for the sector (a) 
reads 



(2^,'='*^'=^) , (x, V) = A(x, 2/) +i?(x, 2/) +C(x, y). 

(a) 

(D29) 



transformed test function have to be analyzed. 

^ a. Pole structure of edge-transformed 
■-test functions In fact, we are interested in 

■^he pole structure of the scaled and then 
translated edge-transformed test functions 
iTxJ^i/sQ^^''^''^f){x',y'), from equation (|83l 



The poles and also the residues of these functions 
are however easily calculated from the poles and 
residues of (Q^^'^^"^ f)ix,y). This is the crucial 
advantage of translational and scale invariance of 

For example, in sector (a), the following poles 
of {Q''^'^'^'^^ f){x,y) with respect to x are obtained: 
xi_2 = ±2i, X3^4 = sy ± 2i, xg^g = -22/ ± 2i, 
2y ± 2i. The corresponding residues are 

1 r - ?— r — ^ ■ -l=Fiy 

27r2jyi '3,4 — 27i-2y' '5.6 — 4^2 (j^2 + i)j^, 

i3ii2L The resulting residues of 



2^7,8 

r7.8 



1 

4Tr2 



{Tx^i/sQ^^'"""f)ix',y') are then given by 



— - ■ subs 



y 



y'-Y 



;niy) 



(D30) 



where the center of mass of the test function is 
X = {X,Y). They are associated to the poles 

of [Tx^A-i/sQ^^'^'^''^ f){x',y'), which are similarly 
given by 



X + S ■ subs 



y 



Y 



(D31) 



In order to obtain the true residues with respect 
to the x'-contraction with the Poisson kernel, one 
only has to evaluate the Poisson kernel at the 



poles ( 48 ) and multiply r[ with the value. 
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b. Pole structure of Poisson kernel The pole 
structure of Vladimirov's Poisson kernel is rather 
straightforward to compute, however rather 
lengthy expressions result for the poles. Similar 
to above, the residues have to be multiplied by 
the function values of the edge-transformed test 

function {Tx^i/sQ^^'^^''''f){x',y') at the pole of 
the Poisson kernel. Poles of the rotated Poisson 
kernel Vr,i}{^m, ^71, x' ,y') with respect to x' are 



■-1,2 



se 
ce - 



—y ±1- 



£112 



se 



V ±1 



se + c 
. m - e»72 



(D32) 
(D33) 



We introduced the short-hand notations 771 = 
(fimC — oJnS, 772 — ifimS + uJnC. The associatcd 
residues of Pr,,? are easily determined. 



c. Residue sum for the x' integral An alge- 
braic expression for the residue sum correspond- 
ing to the x' integral can be generated symboli- 
cally by evaluating 



dx' 7',,4^„„a;„;a;',y') • {rxt^x!,Q^t^''^f){x\v') 



"Poisson" 



"cdec" 



UJn;Xi,y') ■ r'i\ ■ 6'(Ima;-) 

"Poisson" 

+ R427ri.(rxAi/.Qi^'^'=V)(^U')- 



(D34) 



However, it yields rather lengthy formulae, be- 
cause the integrations with respect to x' have to 
be done separately, for each of the sectors ( D23 1 



to ( D28 ) . The growth in complexity is also due to 



amount of parameters which increased dramati- 
cally by introducing the Poisson kernel and in- 
serting the translations 73f and scaling Ai/s of 
the edge functions. It is possible to export the ex- 
pressions resulting from (D34| from the computer 



algebra system to a file of Fortran code which is 
500 kilobytes in size. Similar to the procedure 
described in section [P 1 b| for the second integral, 
the expression ( D34 ) includes removable disconti- 
nuities. 



2. Numerical quadrature of the fourth 
integral 



Due to the vast complexity of expression ( D34 ) 
the last remaining integral 



/4(v3m,w„)= / dy'/3((^,„,w„;y') (D35) 



is evaluated numerically, making use of the 
exported Fortran code. Adaptive integration 
routines from the GNU Scientific Library are 
imposed.'sH Because the resulting matrix elements 
give rise to an inverse problem, it is compulsory to 
achieve a high integration accuracy. By definition, 
the integrand is most distinguishly structured in 
the area y' w F, on a scale 6. Special attention 
has to be drawn to the appropriate integration of 
this range. 

The high-frequency tails (—00, —R] and [R, 00) 
need to be integrated out separately, where R is 
the integration range of the conventional quadra- 
ture. For some parameter values of 6, A, etc., 
problems with the convergence of these high- 
energy integrals may occur, due to floating point 
precision. Choosing a finite interval extending to 
±max(10^ \X\ ■ 10^, |r| • lO^) is then usually suf- 
ficient for numerically satisfactory data. 

A typical shape of the integrand I3 is shown in 
figure 25 The structure at y' « F is not necessar- 
ily 5-shaped, but depending on the values of cpm 
and LOn it may rather look like the Hilbert trans- 
form of such. The integral may be computed at 
each point (A, Y) of the A discretization lattice 
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-0.05 




of the inferred (i? = 
enceilOlfor details. 



0)-edge function. See refer- 



FIG. 25: Example for as a function of y' , us- 
ing (fim = 2.0, LUn — 2.0, with test function loca- 
tion X = 0.0, Y = 4.0, and test function width 
5 = 10"'^. The interacting branch cut geometry is 
used for the determination of r and 



for all values of the simulation data (iipm, ii^n) on 
a computer cluster. In practice, the computation 
of the matrix elements has to be done only once 
for each temperature (3, regardless of the bias volt- 
age. This is because an adjustment of the A grid 
is not necessary in the latter case. 

In future applications, one should aim at sym- 
bolically programming the residue sum of the 
fourth integral and then take the limit (5 — > an- 
alytically. 



Appendix E: MaxEnt implementation for 
data from multiple wedges 

In this appendix chapter, the implementation 
of the MaxEnt algorithm for the Q-mapping is 
described. Details on the computation of the nu- 
merical representation of Q were provided in ap- 
pendix [D| 

The local test function width for the map 
Qr,^ can be adjusted to the local grid resolution 
when the function A{x) is discretized. The inverse 
problem for the inference of spectral properties 



using assumption (65) is, by construction 



lmG{iipm,ii^n) = {Qr^i}A){iipjrt,ii^n)- (El) 

The values r, are those which specify the T'-^'' '' 
branch of G the vector {icpm, is located in, as 

defined by point 1 of section |IIID| The spectral 
function of the dot electrons can then again be 
gained by evaluating along the physical line. 



a. Discretization of A{x^ 



The single-wedge MaxEnt-based analytic con- 
tinuation problem proposed in Ref. \W\ only re- 
quired a rather straightforward discretization of 
the function A{x^, Xi_j). In contrast, for the mul- 
tiwedge mapping, the discretization of the edge 
function has to pay tribute to the strong inter- 
twining of edge structure and branch cut struc- 
ture which is revealed by Q^^'^so) (^f. paragraph 



IVDOh) 



Especially, it turns out that the limiting be- 
haviour along the singular directions of A has to 
be captured numerically. In terms of the multi- 
wedge approach, the singular directions dominate 
the mathematical structure. In our experience, 
also the lateral structure of A{x) along the singu- 
lar directions has to be resolved. We constructed 

(i) (1) 

a grid as follows. Let x\ and be two variables 
which are discretized on i-th and j-th logarithmic 
grid points around zero, respectively. Then the 
VI 8 10^ 
20 I _4 5 



grid ay- 



yields an appro- 



(E2) 



priate discretization of the edge, because the given 
matrix maps the double-cone ]R+ x M+ U x R~ 
and its complement to the wedges defined by the 
singular directions. The numerical test function 
width Sx can then be adjusted to the local grid 
resolution. 

Also the high-energy structure of the Green's 
function has to be taken into account explicitly, 
because along the singular directions it does not 
decay. In practice it seems to be important to 
have a very large logarithmically discretized fit 
region, for which in practice a x^/T region of at 
most [—800, 800] is subject to modifications by the 
MaxEnt algorithm and a x^j/T region of at most 
[—400,400]. The singular-direction contributions 
beyond this range also prove not to be negligible, 
in a test with the bare Green's function (see also 
the Go benchmark below). In order to take them 
into account, their contribution up to very large 
energies {x^ « lO^F) is computed assuming a Go- 
like structure along the directions, positioning ad- 
equately weighted (5-spikes along it and substract- 
ing the corresponding contributions from the raw 
data, as done for the negative-spectral-function 
contributions of static observables in the first pa- 
per. 
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Kernel structure 



The kernel Qr,t5 raay exhibit rather sharp struc- 
tures with respect to the A function space. In 
particular, this may be the case in regions where 
the to-be-determined A is expected to be very 
smooth and physically noninteresting. Conse- 
quently, for these regions, the MaxEnt discretiza- 
tion grid would be chosen rather coarse-grained. 
These potentially disturbing structures can al- 
ready be seen from the formal structure of Vr.^ 
which features strong anisotropics. The con- 
volution with the (transformed) test functions 
Q^edge) j^^^ is in general no cure for this problem, 

because 2^°'^^°'' /x.i even more sharply struc- 
tured, on the scale S, which is of the order of the 



discretiziation scale (see figure 17 1 



In order to discuss this in more detail, some 
matrix elements of Qr.-g are plotted in figure 26 
The orientation of the considered data point in 
Matsubara space defines the orientation of the 
structure which emerges in the kernel with respect 
to the (x^p^x^) coordinates of A-space. A ma- 
jor qualitative difference to the structure of the 
single-wedge kernels Vr.^ is the emergence of dis- 
tinguished negative regions. They are generated 
by the combinations of Hilbert transforms within 
the edge-to-edge map Q'"^'^^^\ As such, they are 
a direct consequence of the branch cuts. The neg- 
ative and positive regions spread over a compa- 
rably wide range and will compete in the process 
of Bayesian inference, in which several (itpm, ia;„) 
pairs and differently overlapping combinations of 
positive/negative regions are involved. The wide 
range of the regions appears to result from super- 
imposing the l/x tails of Q^°'*^°''/x,i5 which are 
dominant for tt/2 and -d w 3tt/2 and absent 
for I? R:! 0, as well as ?9 ~ tt. Note that since 



the continuity assumption ( 65 1 becomes exact for 



larger energies, this feature can be expected to be 
contained in the kernel of an optimal continuation 
theory of Green's functions within the Matsubara 
voltage formalism. 

The kernel structure moreover indicates that 
due to the leverage of the single-wedge con- 
straint (as it applied to the MaxEnt calculations 
in Ref. llOp , the nonequilibrium spectral function 
could now well be resolved. In the following, the 
interacting branch cut geometry will always be 
used for the operator Qr,i?- For brevity, the ac- 
cordingly defined operator will be shortly written 
as Q, since r and d are now well-determined. 



^0 



location within 
wedge geometry j 




-4 4 8 

(a) m = 1, n = (/3r = 5.0) 



^0 




(b) m = 1, n = 3 {I3r = 5.0) 

FIG. 26: (color online) Cut through Qr,i? h^ 
A-space for different pairs of iy^ and ia;„, at 
PV — 5.0. The wedge opening ratios r and 
wedge orientations d (figure 15 1 are chosen ac- 
cording to the interacting branch-cut geometry. 
The "nonequilibrium" line represents the location 
of the dot-electron spectral function for a system 
with source-drain voltage e$ = 0.5r. Wiggly 
structures at higher energies result from the in- 
creasingly coarse-grained A grid. 



c. N on-interacting Green's function as benchmark 

The fundamental assumption of this chapter, 
equation (65), is exact for Go (equation (55)). As 
a consequence, we use the noninteracting Green's 
function as a benchmark for our multiple-wedge 
numerical analytic continuation procedure, al- 
ready assuming the interacting branch cut geome- 
try for the construction of Qr.§, which is certainly 
also valid for Gq. At present, from a numerical 
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FIG. 27: Successive improvement of kernel qual- 
ity by averaging out the local kernel structure 
within the local A grid resolution. Data are shown 
for = 0, $ = 0, ;9r = 5, 71 0,...,9, 
m — —3, . . . , 3 and a realistic mock diagonal co- 
variance matrix C = diag ( ^°p2 ) ■ Abscissa de- 
notes x^/r, ordinate denotes x^/T, grayscale de- 



notes A{x,^ 



in units of F ^ . 



point of view, the method is comprised of two tech- 
nically challenging consecutive steps: First, the 
kernel and its high-energy convolution with the 
Green's function have to be evaluated numerically 
up to a certain precision. Second, an appropriate 
default model has to be defined and the MaxEnt 
must converge to a good estimate in a controlled 
way. 

In order to test the performance of the first step, 
we can take the exact solution as default model 
and run the MaxEnt with the discretized kernel. 
By construction, due to the design of Bryan's 
algorithnP, MaxEnt changes of the A{xip,Xu:) 
function will directly correspond to the numeri- 
cal errors in the computation of the kernel ma- 
trix elements: evidence for changes of the exact 
solution is taken from the exact data due to nu- 
merical imperfections in the kernel. Without in- 
tegrating out the sharp structures of the kernel 
properly for these regions, serious artifacts are ob- 
tained even for larger test function broadnesses Sx- 
This can be seen in the "historic" MaxEnt data 
shown figure 27 Here, is chosen adaptively 
with respect to the local kernel resolution, namely 
Sx = 0.3- (local kernel grid resolution). The Max- 
Ent is able to modify the A function on a large 
grid varying over the ranges x^p G [—800,800], 
Xuj & [—400,400]. As the local kernel resolution is 
increased, averaging out its structure within the A 
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FIG. 28: Sensitivity to the default model, using 



the same parameters and scales as in figure 27d 



grid, an increasingly appropriate discretization of 
the kernel is obtained. In the computations shown 
in figure [27j realistic covariance weights for the 
imaginary-time data were assumed. If numerical 
errors were included into the realization 

of the kernel, there would probably be stronger de- 
viations from Ao{xip, x^) than observed. For some 
single points the Sx is so small that the adaptive 
quadrature of the fourth integral in Q does not 
converge. This can be seen best in figure [27d| be- 
cause here the kernel discretization grid is eight 
times finer than the A{xip, Xi^) discretization grid. 
Similarly, in the nonequilibrium situation, $ 7^ 0, 
the function Aq is not significantly altered by the 
8 X 8-averaging kernel. This was tested explicitly 
also for large bias voltages, such as e<I> = F. 

The performance of the second step can be 
tested by using a kernel realization which suc- 
ceeded in the first step and then performing runs 
with a modified default model. Because the non- 
interacting Aq{x^,x^) function has the correct 
singular behaviour as x — >■ cxd we investigate the 
dependence of the MaxEnt results on the follow- 
ing default models: 



(E3) 

As compared to Ao{x^p, x^), the width of the 
Lorentzians is varied. Using the best-quality ker- 
nel, i.e. 8x 8-averaging (see figure 27d), increasing 
the default-model width quickly results in spuri- 
ous features in the low- to intermediate-energy re- 
gion, even though Q represents an exact relation 
between data and A and the numerical represen- 
tation of Q is sufficiently accurate. Away from 
the low-energy region also for a = l.SF a good 
agreement with Aq is obtained, i.e. a sharpened 
structure along the cross-shaped directions with 
an approximately correct amplitude (as compared 
to figure 27a). The strong sensitivity of especially 



the low-energy range on the default model may 
be interpreted as a result of the subtle interplay 
of positive and negative regions of high- amplitude 
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kernel matrix elements for different (i(/3„i, iw„). 
The structure of the matrix elements was dis- 
cussed above and plotted in figure [26) As shown 
in figure [29] a problem often encountered for not 
well-chosen default models is apparently an in- 
crease of spectral weight in the low-energy region 
1 3; I 0, which exceeds the color scale used in the 
plots by up to a factor of three, even for moderate 
deviations of CTdof from T. This is unfortunate, be- 
cause not only for spectral functions unphysically 
high values may be deduced, but also the over- 
all weight of the spectral function is too large. 
However, since the kernel Q imposes an exact re- 
lation on Go and is resolved well enough, this un- 
fortunate aspect is identified as a pure MaxEnt 
{^^ second step^^) artifact. As such it is no concep- 
tual problem of the Q-approach and can in prin- 
ciple be removed by developing a more sophisti- 
cated MaxEnt algorithm which imposes the phys- 
ical constraints as prior information. In fact, this 
issue can be significantly reduced by a careful but 
straightforward analysis of the posterior probabil- 
ities within a set of smooth default models. 
From our data we can conclude that de- 



fault models with the shape (E3) are appar 



ently not of much use for functions whose high- 
energy behaviour along the singular directions is 
a Lorentzian with width T. Once the high-energy 
structure is known to be such, an interesting ex- 
periment is to flatten out the low- to intermediate- 
energy structure of the default model, by impos- 
ing an x-dependent 



<Jdcf{x^) = r + (o-dcf - r)- 



i?2 



(E4) 



where R is the flattening radius and CTdof ^ T 
is a strong flattening of the default model's low- 
energy region. For Gq it turns out that the re- 
sulting MaxEnt solution is practically identical to 
the (Tdcf = r solution. Consequently, the "second 
siep" -artifact for Go of overshooting low-energy 
values (Fig. 29 ) can just be cured by imposing 
the correct high-energy limit. The low-energy ar- 
tifact is thus caused by missing a-priori informa- 
tion about the high-energy structure. This ap- 
pears to be another manifestation of the fact that 
the kernel Q puts a large range of energy scales 
in relation to each other. 



1. Application to the interacting model 

Switching on a finite Coulomb interaction, one 
has to be aware of the fact that the Q-mapping 
can no longer be expected to be fully exact. How- 
ever, a special case of the assumption, namely the 
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FIG. 29: MaxEnt results for flat low-energy de- 
fault models (E4|, R — 5. As compared to 
Fig. |28b[ the quality of low-energy data is in- 
creased significantly, due to the correct high- 
energy behaviour of the default model. Scales are 



as in figure 27d 



fitting ansatz in Ref. [1, is found to yield reason- 
able results which agree with other methods up 
to a certain exteniPS. Therefore, it seems worth- 
while to investigate how far one can go with the 
controlled MaxEnt approach to the inversion of 
the Q-mappin^^ll 



a. Lateral structure along singular directions 

As shown in the preceding section, the a- 
posteriori determination of a most adequate ap- 
proximate a-priori picture of the high-energy 
structure is crucial for the success of the Max- 
Ent procedure. As sketched in figure [30| at fi- 
nite U, one may, for example, expect the lateral 
structure be an unphysical copy of a spectral func- 
tion, i.e. two Hubbard peaks with possibly an ad- 
ditional peak associated to a quasi-particle res- 
onance. Such a structure would extend over a 
range « U. However, the two parallel Hubbard 
peaks can be expected to approximately have a 
Lorentzian structure of width « F and would gen- 
erate a type of branch cut in the Q-mapping which 
is equivalent to the one in Go. In the strongly cor- 
related regime, Hubbard satellites may be broad- 
ened up to a width of 2r, due to many-body 
correlation^^. 

It is a priori uncertain to which extent either 
of the intuitive pictures in figure [30] is correct. 
However, one of the conceptual strengths of the 
Q-mapping is the precise rendering of the high- 
energy structure of the imaginary-voltage theory 
(cf. section 
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One can expect that only a 
characteristic width of the lateral structure along 
the singular directions is needed in order to model 
the correct high-energy contribution to the am- 
plitude of the discontinuity of G(z^,Zt^) at the 
low-to-intermediate energy portions of the branch 
cuts. Based on this, we can investigate the poste- 
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rior probability Pr(crdcf |G) for default models (E3 ) 
as a function of their width ddcf- 



AAA AAA A 

r 

(a) nontrivial physics at(b) 
±x^/2 



"asymptotic free- 
dom" 



FIG. 30: At high energies, one might (a) ex- 
pect the lateral structure of A{x^^ x,^) to be com- 
posed of two Hubbard peaks and possibly a quasi- 
particle resonance which combine to the physical 
spectrum A{uj) at the intersection point. In the 
complementary scenario (b) , the function A would 
not differ from the noninteracting one at high en- 
ergies. 




FIG. 31: (color online) Posterior probability of 
the default model (|E3]) at /3r = 5, e$ = T for 
several interactions strengths. The result is found 
to be essentially independent of the bias voltage. 



The thus determined most probable ddof then 
serves as an effective description of the high- 
energy structure A{x^p,Xi^) for the actual com- 
putations. However, as input data from the 
QMC simulations, only low- to intermediate- 
energy data are available. Therefore, the posterior 
probability probe with respect to default models 
(E3) has to be interpreted with care. In figure 



The kernel validated in figure 27d has been used 



31 posterior probabilities for different interaction 
strengths are displayed. Due to the width being 
significantly larger than 2r for [/ = ISF it is ob- 
vious that the lateral width cannot solely be in- 
terpreted as a signature of the Hubbard bands. 
Merely, the overall Lorentzian broadness of the 
spectral function seems to be obtained. Based 
on our data, neither of the scenarios of figure [30] 
can be preferred. However, based on our expe- 
rience, the most probable high-energy structure 
also yields reasonable results in the case of com- 
parably strong interactions. Thus, in the practi- 
cal computations, first the most probable default 
model is identified. As a next step, the actual 
spectral functions are estimated. 
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Sheaves are very general mathematical concepts, 
but here one may just imagine the set of holomor- 
phic functions f/ — > C carried by a given domain 
U C C*. 
See section 

Because in our application e will not necessarily be 
small, we changed to the notation r. 



For the definitions see figure 15 



To be more precise, for technical reasons (conser- 
vation of the spectral function normalization for 
any causal selfenergy) the fits in Refs. 01251 were 
performed with respect to the selfenergy, not with 
respect to the Green's function. Because the Q- 
mapping is also exact for ansatz functions of the 
fits presented in Ref. [T] the fits are essentially spe- 
cial cases of the present work, though. The first- 
order Pade approximant fits presented in Ref. [23] 
represent a different kind of generalization. 



